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ABSTRACT 

some  applications,  it  is  desirable  to  have  an  antenna  which  can 
adaptively  form  antenna  pattern  *hulls*  in  the  directions  of  a  collection  of 
interference  sources.  These  nulls  are  generally  required  to  be  below  some 
specified  level  of  directive  gain  over  the  entire  system  bandwidth.  Typically 
the  nulls  are  formed  by  adaptively  weighting  each  output  of  an  N-port  antenna, 
such  as  an  N-element  array.  For  large  arrays  operating  over  wide  percentage 
bandwidths,  some  sort  of  frequency  compensation  is  required  at  each  output 
port  to  accomplish  broadband  nulling.  One  technique  commonly  used  is  to  employ 
a  •'tapped*  delay  line  at  each  element  output,  with  controllable  weights 
(frequency  independent)  at  the  output  of  each  tap.  The  objective  of  this  report 
is  to  develop  some  insight  into  the  way  in  which  delay  line  compensation  leads 
to  improved  performance,  to  develop  some  quantitative  estimates  of  how  perfor¬ 
mance  varies  with  antenna  and  delay  line  parameters,  and  to  develop  some  tools 
useful  in  the  performance  evaluation  of  frequency  dependent  weighting. 

Practical  considerations  dictate  that  only  a  minimum  number  of  taps / 
element  be  utilized.  A  primary  consideration  is  to  minimize  the  number  of  taps 
required  for  a  specific  system  performance.  With  this  objective  in  mind,  we 
develop  an  analytic  solution  characterizing  the  two  element  array,  employing 
a  two-tap  delay  line  with  each  element.  This  provides  valuable  insight  into 
the  mechanism  of  delay  line  compensation.  ^Results  for  this  simple  case  will 
then  be  extended  to  the  N-element  planar  array  using  Nt  taps/element.  That 
tap  spacing  which  optimizes  the  array  performance  from  the  viewpoint  of 
providing  wideband  nulls  on  sources  incident  on  the  array  from  locations 
arbitrarily  specified  over  the  antenna  field  of  view  is  developed,  and  applied 
to  numerous  types  of  array  configurations.  The  composite  performance  of  nine 
specific  array  configurations  is  used  to  estimate  the  number  of  taps  and  tap 
spacing  which  would  be  required  for  an  arbitrary  array  configuration,  as  a 
function  of  the  array  and  delay  line  parameters  Ay^  =  2w  ®  sinB  ,  FBW, 

and  Nt,  where  FBW  is  the  fractional  bandwidth  (100*FBW  =  percentage  nulling 
bandwidth)  0^  is  the  maximum  scan  angle  off  broadside,  4>q  the  tap  spacing 
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(in  degrees  relative  to  the  center-band  frequency)  D  is  the  maximum  aperture 
dimension,  and  N£  the  number  of  taps.  The  implications  of  using  an  Applebaum- 
Ho wells  type  of  adaptive  processor  to  determine  the  delay-line  output  weighting, 
and  the  algorithm  effects  on  performance,  are  also  developed. 
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Fig.  8.  Cancellation  (dB)  vs  Ps^j^  for  the  two-element,  two  tap 
array  with  Ay  •  FBW  and  tap  spacing  fixed:  a)  Ay  •  FBW  =1.0 
b)  Ay  «  FBW  =  4.0. 
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I. 


INTRODUCTION 


For  some  applications,  it  is  desirable  to  have  an  antenna  which  can  adap¬ 
tively  form  antenna  pattern  minima,  or  "nulls"  in  the  direction  of  a  collection 
of  interference  sources.  These  minima  are  generally  required  to  be  below  some 
specified  level  of  directive  gain  over  the  entire  system  bandwidth.  When  this 
bandwidth  is  quite  large  (for  example,  when  bandspreading  is  used  in  the  modu¬ 
lation  format  of  a  communications  signal),  the  level  of  nulling  which  can  be 
achieved  becomes  dependent  on  the  aperture  size  (in  wavelengths),  the  antenna 
field  of  view  (FOV) ,  and  the  RF  fractional  bandwidth  (FBW) .  Typically,  such 
nulls  are  formed  by  adaptively  weighting  each  of  the  outputs  of  an  N-port 
antenna  such  as  an  N-element  array,  or  N-beam  multibeam  antenna,  and  then 
summing  these  weighted  outputs  so  as  to  form  minima  in  the  radiation  pattern  in 
the  directions  of  the  interference  sources,  and  relative  maxima  in  the  directions 
of  the  desired  users  of  the  system.  For  large  antennas  operating  over  wide 
percentage  bandwidths,  some  sort  of  frequency  compensation  is  required  at  each 
antenna  output  port  to  accomplish  broadband  nulling.  One  technique  commonly 
used  is  a  frequency  dependent  weight  consisting  of  a  "tapped"  delay  line  with 
controllable  weights  (frequency  independent)  at  the  output  of  each  tap,  and 
the  weighted  tap  outputs  summed.  The  net  effect  is  to  adaptively  synthesize 

the  "optimum"  frequency-dependent  weight  at  each  element.  Several  papers  have 

(1  2  3) 

addressed  this  problem  .  The  objective  of  this  report  is  to  develop  some 

additional  insight  into  the  way  in  which  delay  line  compensation  leads  to 
improved  performance,  to  develop  9ome  quantitative  estimates  of  how  performance 
varies  with  the  antenna  and  delay  line  parameters,  and  to  develop  some  tools 
useful  in  the  performance  evaluation  of  frequency  dependent  weighting. 

Practical  considerations  dictate  that  only  the  minimum  number  of  taps  per 
element  required  to  achieve  the  desired  performance  be  utilized.  For  an  N- 
element  array  using  taps  (Nt=l  corresponds  to  frequency  independent  weights 
with  no  delay  lines),  then  N*Nt  controllable  weights  are  required.  A  primary 
consideration  is  to  minimize  required  for  a  specific  system  performance. 

With  this  goal  in  mind,  it  is  useful  to  consider  the  performance  characteristics 
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of  an  N-element  array  using  frequency  independent  weighting,  so  that  the  need 
for  incorporating  the  more  complex  delay  line  processor  can  be  quantified. 

Fig.  1  illustrates  an  estimate  of  the  cancellation  performance,  using  frequency 
independent  weighting,  for  an  N-element  linear  array  for  the  cases  N=2,  7  and 
17.  Each  array  case  forms  N-l  nulls  on  N-l  interference  sources  separated 
greater  than  a  half-power  beamwidth  (HPBW)  from  its  neighbor.  For  our  puposes, 
the  cancellation  C  plotted  in  Fig.  1  is  defined  to  be  the  quotient  of  the 
interference  to  thermal  noise  ratio  at  the  array  output  after  nulling  to  that 
before  nulling  (see  Eqs.  38  and  39  of  the  text),  when  the  interference  is 
assumed  to  consist  of  white  noise  over  the  bandwidth  BW.  The  quiescent  radia¬ 
tion  pattern  of  the  array  is  assumed  to  be  "earth  coverage"  (i.e.,  a  single 
array  element  is  used  as  the  reference  element,  where  the  element  pattern 
covers  the  earth  FOV  with  applications  pertaining  to  an  uplink  to  a  satellite 
at  geosynchronous  orbit).  For  convenience  we  define  the  fractional  bandwidth, 
FBW,  as  BW/fp,  where  f^  is  the  operating  frequency  and  BW  is  the  nulling  band¬ 
width  (in  hertz).  The  development  leading  up  to  Fig.  1  is  presented  in  Appen¬ 
dix  A.  Referring  to  the  figure,  it  is  convenient  to  characterize  the  array 
performance  according  to  the  variables  D/X,  sin©m  and  FBW,  where  D/X  denotes 
the  aperture  diameter  in  wavelengths,  and  ©m  is  the  maximum  scan  angle  defining 
the  field  of  view  over  which  the  array  must  operate.  The  cancellation  degrades 
rapidly  as  any  of  the  parameters  D/X,  ©m  or  FBW  is  increased.  More  precisely 

C  a  (Ay  •  FBW)2,  where  Ay  =  2tt  •  •  sin©  .  If  the  desired  cancellation 

m  m  A  in 

performance  is  specified.  Fig.  1  can  be  used  to  estimate  the  corresponding 

maximum  allowable  value  of  Ay  consistent  with  this  performance  specification. 

m 

If  system  requirements  dictate  that  a  value  of  Aym  greater  than  this  maximum  is 
required  (for  example,  the  minimum  allowable  user-interference  separation 
dictates  the  choice  of  D/X),  frequency  dependent  weighting  requiring  Nt>l 
must  be  used.  As  an  example  of  the  use  of  Fig.  1,  if  we  desire  C  5  -30  dB, 
then  Ay^'FBW  must  be  less  than  0.25  if  frequency  independent  weights  are  used. 
For  the  geosynchronous  satellite  uplink  geometry,  0^  is  fixed  at  approximately 
9°.  D/X  is  determined  by  the  user- interference  resolution  required.  The 
tolerable  bandwidth  with  frequency  Independent  weighting  is  then  specified  by 
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the  requirement  AYm*FBW  ^0.25.  Results  similar  to  Fig.  1  specifying  the 
performance  of  the  delay  line  processor  for  use  with  specific  array  configu¬ 
rations  will  be  presented  in  Sections  III  and  IV. 

Observe  from  Fig.  1  that  the  cancellation  estimates  are  relatively  inde¬ 
pendent  of  the  number  of  array  elements  used.  It  will  be  seen  later  that,  in 
fact,  the  estimates  of  Fig.  1  serve  as  an  upper  bound  on  the  cancellation 

performance  of  the  more  complex  planar  arrays.  However,  the  general  behavior 

2 

indicating  C  a  (Ay^  •  FBW)  remains  valid.  Only  the  constant  of  proportion¬ 
ality  changes  with  array  configuration.  It  should  be  pointed  out,  however, 
that  this  dependence  is  characteristic  only  of  array  antennas.  For  multiple 

beam  antennas  (MBA's)  employing  beams  formed  using  time-delayed  techniques,  it 

2 

can  be  shown  that  C  a  FBW  ,  and  is  relatively  independent  of  D/A  and  0  .  This 

m 

results  from  the  time-delayed  nature  of  each  of  the  beams.  A  detailed  discus¬ 
sion  of  the  MBA  bandwidth  behavior  is  presented  in  Ref.  4.  Because  of  these 
basic  differences  between  the  MBA  and  the  array  antenna,  the  results  charac¬ 
terizing  the  array  configurations  treated  in  the  following  sections  cannot  be 
directly  extrapolated  to  the  MBA. 

In  the  following  section,  we  will  introduce  several  concepts  which  are 
fundamental  to  obtaining  a  quantitative  assessment  of  the  performance  of 
adaptive  antennas  using  a  delay  line  processor.  Direct  assessment  of  the  band 
with  performance  of  an  N-element  antenna  used  in  conjunction  with  an  adaptive 
processor  is  difficult  unless  the  array  is  required  to  form  N-l  independent 
nulls.  This  is  because  the  array  possesses  N-l  available  degrees  of  freedom, 
and  forming  less  than  N-l  nulls  allows  the  processor  to  compensate  for  any 
narrowband  tendency  inherent  in  the  antenna  by  using  the  additional  available 
degrees  of  freedom  to  form  wide-band  nulls.  However,  this  type  of  wide-band 
null  tends  to  be  spatially  broader  than  the  corresponding  narrowband  case,  so 
that  the  spatial  resolution  of  the  array  is  compromised.  When  a  tapped  delay 
1  oe  is  used  to  wide-band  the  array,  the  narrowband  spatial  resolution  of  the 
array  is  preserved.  To  avoid  this  ambiguity,  and  to  establish  a  basis  for 
comparing  array  performance  with  and  without  a  delay  line  processor,  we  intro¬ 
duce  the  concept  of  "orthogonal"  sources.  By  definition,  such  sources  use  up 


complete  degree  of  freedom  per  source  to  be  nulled.  Subjecting  the  array  to 
N-l  orthogonal  Interference  sources  necessarily  requires  all  the  available  N-l 
degrees  of  freedom  of  the  array  to  be  allocated  to  spatial  nulling,  with  none 
remaining  for  bandwidth  compensation.  The  improvement  in  performance  realized 
using  delay  line  processing  compared  to  frequency  independent  weighting  can 
then  be  quantitatively  assessed  for  this  type  of  scenario. 

A  second  useful  concept  offering  insight  into  how  delay  line  compensation 
leads  to  improved  performance  is  that  of  "optimum  weighting".  The  set  of 
optimum  weights,  by  definition,  has  the  precise  frequency  variation  required  to 
perfectly  null  all  the  N-l  orthogonal  interference  sources  over  the  entire 
bandwidth  of  interest.  The  frequency  dependence  of  these  optimum  weights  can 
readily  be  determined.  The  function  of  the  delay  line  processor  is  to  synthe¬ 
size  an  effective  weight  variation  with  frequency  which  approximates  this 
optimum  weight  variation.  The  degree  to  which  this  is  possible  depends  on  the 
delay  line  parameters  i.e.,  the  tap  spacing  and  the  number  of  taps.  A  novel 
graphical  presentation  illustrating  this  synthesis  of  the  optimum  weights  via 
delay  lines  will  be  developed  in  Section  II  and  used  throughout  in  evaluating 
the  characteristics  of  specific  antenna  configurations. 

The  performance  of  an  N-element  adapative  array  is  strongly  dependent  on 
the  control  algorithm  used  in  the  delay  line  processor.  In  Section  III  we 
present  an  analysis  of  this  dependence  assuming  that  an  Applebaum-Howells  type 
adaptive  feedback  processor ^  controls  the  weights  at  the  output  of  each  tap 
of  the  delay  line.  For  this  type  of  processor,  the  quiescent  radiation  pattern 
of  the  array  is  governed  by  "beam-steering"  controls  applied  to  the  control 
loops.  Upon  adaption,  the  weights  are  set  to  a  value  which  reduces  the  radia¬ 
tion  pattern  in  the  direction  of  the  interference  sources  to  a  minimum  level 
dependent  on  the  loop  gain  (or  processor  dynamic  range) ,  while  simultaneously 
yielding  the  best  rms  approximation  to  the  quiescent  radiation  pattern  in 
directions  away  from  the  interference.  Choice  of  loop  gain  (i.e.,  the  effec¬ 
tive  power  level  of  the  interference  sensed  by  the  processor)  determines  the 
tradeoff  between  depth  of  null  and  rms  pattern  error  of  the  adapted  pattern 
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relative  to  the  quiescent  pattern.  The  dependence  of  the  array  performance  on 
processor  parameters  will  be  evaluated  for  the  two-element  array  using  two 
taps/element.  For  this  simple  case,  detailed  analytical  solutions  can  be 
obtained,  providing  a  physical  feeling  for  how  delay  line  compensation  improves 
the  array  nulling  bandwidth.  Proper  choice  of  the  steering  vector  required  to 
keep  the  array  from  "turning  itself  off"  will  also  be  developed.  It  will  be 
shown  that  only  in  the  limit  of  very  large  loop  gain  is  the  full  potential  of 
the  delay  line  processor  realized. 

Finally,  in  Section  IV,  the  performance  of  specific  N-element  linear  and 
planar  arrays  operating  over  a  fixed  FOV  will  be  quantified  as  a  function  of 
Ay^,  the  tap  spacing,  number  of  taps  and  the  nulling  bandwidth  for  processors 
having  very  large  loop  gain.  In  this  limit,  array  performance  becomes  indepen¬ 
dent  of  the  control  algorithm,  and  only  dependent  on  the  array  and  delay  line 
parameters.  Of  prime  Interest  will  be  to  characterize  how  the  cancellation 
over  FBW  improves  as  the  number  of  taps  increases,  so  that  only  the  minimum 
number  of  taps  might  be  used  for  specific  system  applications.  The  composite 
results  determined  from  numerous  array  configurations  will  be  used  to  estimate 
a  specific  value  of  Nt  required  to  obtain  a  specified  cancellation  performance 
for  an  arbitrary  array  configuration.  These  results  are  summarized  in 
Figs.  24a,  b  and  c  of  the  text. 
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II.  PERFORMANCE  EVALUATION  CONCEPTS 


A.  The  Concept  of  Orthogonal  Interference  Sources 

As  discussed  in  the  previous  section,  it  is  useful  to  introduce  the  con¬ 
cept  of  orthogonal  interference  sources.  By  definition,  such  sources  use  up  a 
complete  degree  of  freedom  available  to  the  array  per  source  to  be  nulled. 
Subjecting  an  N-element  array  to  N-l  orthogonal  sources  necessarily  requires 
all  the  available  N-l  degrees  of  freedom  of  the  array  to  be  allocated  to  spa¬ 
tial  nulling,  with  none  remaining  for  bandwidth  compensation.  Mathematically, 
this  condition  can  be  stated  as  follows: 

Define  the  source  direction  vector  jj,  whose  components  are  the  outputs 
of  each  element  of  the  array  when  a  narrowband  signal  of  unit  amplitude  from 
direction  0,  4>  is  incident  on  the  aperture,  according  to 

+jk.  y  sin0(x,  cos<J>  +  y,  sin<)>) 

(u)k  =  e  K  k=l . N  (1) 

where  kQ  =  2t\)\  -  ui^/c,  *  is  the  wavelength,  uiQ  the  operating  frequency, 

c  *  3xl08  m/sec,  D  the  diameter  of  the  circle  circumscribing  the  array,  0,  <J> 

the  angular  coordinates  specifying  the  source  direction  relative  to  array 

broadside,  and  the  x,  ,  y,  are  the  normalized  (relative  to  D/2)  element  place- 

■  2  f 

ment  coordinates  (x^  +  y^  <  1).  Hence,  when  u  is  applied  to  the  array 
elements  as  a  weight  vector,  then  the  array  output  response  is  maximum.  (Note: 
"+"  denotes  "complex  conjugate  transpose.")  A  set  of  orthogonal  sources  at 
positions  u^,  u^  satisfies  the  orthogonality  conditions 

Hi+  •  Uj  =  N6t  i,j  =  1 . N  (2) 

where  6.  denotes  the  Kronecker  delta  function  (6  «  0,  i  4  j;  6.  ,  =  1, 

**J  1*j  i*J 

i  =  j).  Physically  this  occurs  when  the  sources  are  sufficiently  far  apart  so 

that  a  separate  null  (other  than  a  grating  lobe  null)  is  required  for  each. 

At  least  one  set  of  sources  satisfying  (2)  can  be  found  for  most  arrays  of 
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interest,  and  a  set  which  at  least  approximately  satisfies  (2)  can  be  found 
for  all  arrays. 

Consider  then  J  narrowband  incoherent  interference  sources  emanating  from 
the  directions  0^ ,  ,  j«l,  J.  The  correlation  matrix  characterizing  the 

array  element  voltage  outputs,  E^,  is  given  by 

R.  =  E  <  E,  E  *  l  =R,G  u,  u.+  +  ...  +  8.  G  u.  uT+  +  I  (3) 

=A  I  It  q  j  1  e  -1  -1  J  e  -J  -J 

til 

where  Ge  is  the  element  gain,  R^  is  the  interference  power  level  of  the  j — 
source  (relative  to  thermal  noise  referred  to  the  element  input)  incident  on 
the  array  aperture,  E  {  •  }  denotes  the  statistical  expected  value  over  the 
interference  noise  processes  and  I  is  the  identity  matrix.  The  N-l  degrees  of 
freedom  of  the  array  are  characterized  by  the  N 'eigenvalues,  s^,  ...,  sN,  of 
the  N  x  N  matrix  gA.  For  J  sources,  J  of  the  N  eigenvalues  of  gA  will  be 
greater  than  unity,  while  the  remaining  N-J  will  be  1.  However,  if  several  of 
the  interference  sources  are  not  orthogonal,  i.e.,  they  are  in  close  proximity 
(relative  to  the  array  half-power  beamwidth,  HPBW) ,  then  not  all  of  the  s^, 

...,  Sj  will  be  significantly  greater  than  1.  Otherwise  said,  the  array  in 
this  case  is  not  using  a  single  degree  of  freedom  per  interference  source 
present  in  order  to  accomplish  the  nulling.  If  the  sources  are  orthogonal  as 
defined  in  (2),  however,  all  of  the  s^,  ...»  s^  must  be  significant,  and  in 
fact  this  condition  leads  to  eigenvalues  s^  =  R^,  ...,  s^  *=  R^.  Equation  (2) 
defines  the  conditions  for  orthogonal  sources,  and  the  presence  of  N-l  such 
sources  necessarily  uses  up  all  of  the  available  N-l  degrees  of  freedom  of  the 
array.  (We  note  that  if  the  array  element  spacing  is  so  large  that  grating 
lobes  appear  over  the  FOV,  then  the  number  of  possible  positions  for  the  u^  to 
satisfy  Eq.  (2)  is  multiplied.) 

The  concept  of  orthogonal  sources  can  be  quite  useful  in  characterizing 
the  performance  of  an  array  when  operating  over  a  wide  bandwidth.  Consider  the 
presence  of  J  orthogonal  sources  operating  over  a  fractional  bandwidth  FBW. 

The  correlation  matrix  for  incoherent  white  noise  sources  is  then  given  by 
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(4) 


h  =  Ge  R1  <  ^1  % 


+  G  R, 


<Hj  Uj 


+  I 


where  the  bracket  <  *  >  denotes  the  integral 

l+FBW/2 


l-FBW/2 


(5) 


and  w  =  denotes  the  normalized  frequency.  Because  of  this  frequency 

average,  J  sources  now  lead  to  the  presence  of  an  additional  group  of  eigen¬ 
values  Sj+1>  Sj+2«  •••  sj+l  arise  due  to  bandwidth  (see  Ref.  6  for  a 

more  detailed  discussion  of  this  phenomenon).  If  FBW  is  large,  the  sj+^> 

. . .  Sj+L  become  significant  and  forming  a  wideband  null  on  the  sources  uses  up 
additional  array  degrees  of  freedom.  This  has  two  possible  effects:  if 
J+L  <  N,  then  the  array  will  form  broadband  nulls,  but  these  will  also  tend  to 
be  spatially  broader  than  narrowband  nulls  and  so  will  yield  reduced  antenna 
gain  to  the  desired  users  in  close  proximity  to  the  interference  sources;  if 
J+L  >  N,  the  array  cannot  null  all  the  sources.  Both  effects  are  to  be  avoided, 
if  possible,  as  they  result  in  non-optimum  array  performance.  The  case 
J+L  <  N  is  difficult  to  evaluate  quantitatively,  since  the  resulting  loss  in 
gain  to  a  user  signal  which  is  tolerable  is  dependent  on  the  system  link 
margin  parameters  and  the  locations  of  the  user  and  interference  sources. 
However,  if  we  consider  the  scenario  consisting  of  N-l  equal,  orthogonal  inter¬ 
ference  sources,  then  SN_^>:> SN  and  the  resultant  eigenvalue  is  directly 
indicative  of  the  array  bandwidth  characteristics.  Hence  direct  examination 
of  s^,  treated  as  a  function  of  FBW,  offers  a  measure  of  the  nulling  bandwidth 
potential  of  a  particular  antenna  configuration. 


B.  The  Tapped  Delay  Line  and  the  Applebaum-Howells  Processor 

Fig.  2  illustrates  the  general  topology  of  a  tapped  delay  line,  employing 
Nt  taps/element,  used  in  conjunction  with  an  Applebaum-Howells  type  adaptive 
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k  —  ARRAY  ELEMENT 


{ 

WEIGHT  CONTROL 


Fig.  2.  Topology  of gapped  delay  line  weighting  scheme 
illustrated  for  the  k—  element  of  an  N-element  array,  and  using 
and  Applebaum-Howells  type  adaptive  processor. 


processor.  For  simplicity  only  the  schematic  pertaining  to  the  k —  element 
is  shown  in  the  figure.  We  use  the  following  notation:  the  subscript  "k" 
denotes  the  array  element  under  consideration  and  the  subscript  "n"  the  loca¬ 
tion  along  the  delay  line.  4>q  denotes  the  tap  spacing  in  degrees  -  i.e., 

=■  a)  T,  where  T  is  the  time  delay  between  taps.  The  frequency  response 
of  the  kr^-  array  element  is  denoted  as  E^(w)  and  the  frequency  response  of  the 
output  of  the  n^  tap  in  the  k—  channel  as  E.  (w) .  Since  considerable  loss 
may  be  encountered  in  the  power  division  to  the  taps,  it  is  desirable  to 

amplify  before  the  delay  line  input.  This  is  indicated  in  Fig.  2,  depicting  an 

2 

amplifier  at  the  output  of  each  array  element  having  thermal  noise  Hq  .  As  a 
consequence  of  amplifying  before  the  power  split,  the  resultant  thermal  noise 
at  the  output  of  each  tap  in  a  given  channel  is  correlated,  as  contrasted  to 
the  more  familiar  situation,  without  frequency  compensation,  in  which  the 
thermal  noise  is  uncorrelated  between  channels.  Because  of  this  noise  corre¬ 
lation  effect,  and  because  the  algorithm  chosen  is  one  that  tends  to  minimize 
received  power,  it  is  possible  for  the  adaptive  array  to  "turn  itself  off"  if 
adequate  constraints  are  not  Imposed  on  the  weight  control  algorithm.  For  the 
Applebaum-Howells  algorithm,  these  constraints  are  Imposed  via  the  beam-steering 
voltage  V,  Indicated  in  Fig.  2,  which  controls  the  weight  in  the  absence  of 
interference  sources  —  i.e.,  V  specifies  the  quiescent  radiation  pattern  of 
the  array.  The  proper  choice  of  V  will  be  developed  in  the  next  sections,  and 
will  be  seen  to  have  a  substantial  effect  on  the  nulling  performance. 

The  outputs  of  each  of  the  N  •  terminals  of  the  array  define  an  N  • 
x  N  •  N  correlation  matrix  R  which  governs  the  performance  of  the  array.  Once 
R  is  determined,  the  weights  can  be  set  according  to  the  chosen  adaptive 
algorithm^ . 

A  -  [I  +  pR]_1  •  V  (6) 

where  p  is  the  effective  processor  gain,  p  is  a  variable,  dependent  on  hard¬ 
ware  considerations,  and  significantly  affects  the  processor  performance. 
Implicit  in  Eq.  (6)  is  the  assumption  that  N  *  beam  steering  voltages 
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characterized  by  the  vector  V  are  applied  to  each  weight  at  the  output  of  each 
tap,  and  that  N  •  N  adaptively  controllable  complex  weights  are  used.  To 
develop  the  expression  for  R,  consider  the  output  of  the  n^  tap  behind  the 
k— •  element: 


.  -j  (n-lH-w 

^  n  ~  Mw>  e 


(7) 


where  E^Cw)  is  the  array  element  output  response.  The  contribution  to  E^  ^(w) 
from  thermal  noise  results  In  a  noise  covariance  matrix  given  by 


5n 


<  e 


-j  (n-mHgW 


(8) 


The  overall  correlation  matrix  R  can  be  written  aa 

.  *  -j  (n-m)<j>_w 

R  =  n—  <  E  {  E.  (w)  E  (w)  e  }>  +  R 

=  N  K  q 


(9) 


For  J  uncorrelated  white-noise  interference  sources,  and  for  an  array  antenna 
configuration,  this  expression  for  R  reduces  to 


J 

R  =2  Rj  +  Rn  (10) 

j=i 


where 


R1  Ge  1  1+F/>W/2  +J(Yk  j-Y  ,)w  -j(n-m)4>  w 

Sj  (k,q,n,m)  -  -Ijf  ^  J  *>  k’J  q‘j  -  -  (11) 


FBW 


dw 


l-FBW/2 
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0 

where  =  — jj—  sin0j<x^  cosi^  +  y^  ein$j)  denotes  the  center  frequency  phase 

shift  at  the  kth  element  due  to  the  jth  source.  If  R.  is  normalized  relative  to 

2  ^ 

Hq  ,  then  the  resultant  expression  for  R  can  be  written  as 


J 

R(k,q,n,m)  =  ^  SjCk.q.n.m)  +  < 


-j  (n-m)  4>qW 
e  > 


(12) 


This  general  expression  for  R  will  be  used  throughout  the  text  to  develop  the 
specific  properties  of  the  adaptive  array  with  delay  line  processing. 

A  careful  examination  of  the  expression  for  R  in  Eq.  (12)  and  the  adapted 
solution  for  the  weights,  Eq.  (6),  will  indicate  several  parameters  affecting 
the  adaptive  array  performance.  The  most  obvious  are  the  tap  spacing  and  number 
of  taps  per  element.  In  order  to  simplify  a  tradeoff  study  characterizing  the 
''optimum"  tap  spacing,  we  have  assumed  a  uniform  spacing  between  taps.  Clearly 
for  some  scenarios,  this  is  a  non-optimum  choice.  For  example,  for  a  3-element 
array  with  two  sources  Incident  at  specified  angles  of  incidence  it  can  be 
shown  that  non-uniform  tap-spacing  gives  the  best  performance  (see  Section  IV) . 
However,  it  might  be  argued  that  a  fixed  set  of  delay  lines  must  be  designed 
to  handle  arbitrary  angles  of  incidence  over  the  antenna  FOV,  in  which  case 
some  angles  of  incidence  necessarily  yield  better  cancellation  performance 
than  others.  Thus,  optimum  performance  should  be  defined  relative  to,  perhaps, 
the  average  cancellation  performance  over  this  range  of  angles.  In  this  case, 
it  would  appear  that  the  assumption  of  uniform  spacing  between  taps  is  not  so 
restrictive. 

The  algorithm  defined  by  Eq.  (6)  defining  the  adapted  weights  is  power 
level  dependent  as  a  consequence  of  the  effective  processor  gain  p.  Define 
e^  and  a^  to  be  the  eigenvectors  and  eigenvalues  of  R,  respectively.  Then 
Eq.  (6)  can  be  expressed  in  the  form 
N*Nt 

*  ■  £  ‘i/1  %  <13> 

i-1  1 
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Eq.  (13)  specifying  the  adapted  weight  vector  A  is  similar  in  form  to  that  for 
the  more  conventional  Applebaum-Howells  processor  without  tapped  delay  lines. 
Note  that  because  of  the  term  l/(l+us^),  the  eigenvectors  corresponding  to  the 
smaller  eigenvalues  are  weighted  more  heavily,  and  must  contain  the  information 
required  to  null  the  interference  sources.  The  detailed  characteristics  of 
these  eigenvectors  corresponding  to  the  lower  eigenvalues  will  be  developed  in 
the  next  section.  However,  it  is  clear  from  Eq.  (13)  that  the  relative  contri¬ 
bution  of  the  eigenvectors  corresponding  to  the  smaller  eigenvalues  to  the 
adapted  solution  for  A  is  dependent  on  the  processor  loop  gain  y;  i.e.,  only 

eigenvalues  s,  such  that  ys.  ^1  are  sensed  by  the  processor.  Define  sM<v  to 
1  1  MAX 

be  the  maximum  eigenvalue  of  R.  Then  it  can  be  shown  that  s.,.„  *  N  G  R.  for 

=  MAX  e  j 

a  single  source  of  relative  power  level  R^  incident  on  the  array.  Consequently 
|JSMAX  ^etennines  t*ie  dynamic  range  of  incident  power  levels  which  will  be  sensed 
by  the  processor.  Hence  the  performance  of  the  delay  line  processor  and  the 
choice  of  "optimum"  delay  line  parameters  can  be  expected  to  be  power  level 
dependent.  Only  in  the  limit  y SMA^ >  >  1  can  we  determine  optimum  delay  line 
parameters  independent  of  Incident  power  levels. 

C.  Optimum  Weighting  vs  Delay  Line  Synthesized  Weighting 

In  order  to  gain  some  insight  as  to  how  frequency  compensation  obtained 
using  tapped  delay  lines  Improves  the  array  performance,  it  is  useful  to 
define  what  we  might  refer  to  as  an  optimum  weight  variation  with  frequency, 
and  then  examine  how  the  frequency  variation  of  the  delay  line  synthesized 
weight  compares  to  this  optimum.  As  illustrated  in  Fig.  2,  E^(w)  denotes 
the  frequency  response  transfer  function  of  the  k—  antenna  port,  before  delay 
line  processing.  Let  us  define  the  OPTIMUM  weighting,  Aq^(w),  k*l,  ...,  N, 
such  that  when  A^^(w)  is  applied  as  weight  to  the  k—  antenna  port  output,  then 
the  interference  sources  are  ideally  nulled  over  the  entire  band  defined  by 
FBW.  The  output  Vq(w)  as  a  function  of  w  is  given  by 
N 

VQ(w)  ”  2  Aok*(w)  Ek(w) 
k«l 
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(14) 


B 


Noting  Eq.  (14),  the  optimum  weights  Agk(w)  can  readily  be  determined  as 
follows:  define  the  "narrowband"  correlation  matrix,  R^w),  to  be  the  N  x  N 
matrix  at  the  array  element  outputs  when  J  incoherent  narrowband  interference 
sources  at  frequency  w  are  Incident  on  the  array,  and  w  is  considered  a  para¬ 
meter  over  the  interval  1  -  FBW/2  <  w  <  1  +  FBW/2.  Then  Aq^(w)  are  given  by 

Aq(w)  =  [I  +  wR^w)]'1  •  V,,  (15) 


where  V„ 
w,  Aq(w) 
by  FBW. 
using  Nt 
VQ(w)  is 


is  an  Nxl  equivalent  steering  vector.  When  considered  as  a  function  of 
minimizes  the  array  output  power  over  the  entire  nulling  band  defined 
For  comparison,  consider  the  delay  line  synthesized  weights  obtained 
taps/element  for  an  N-element  array.  From  Fig.  2,  the  output  voltage 


given  by 


N 


V">  -  2 

k-1 


__  * 

2  Vn 


n=l 


-j (n-l)$0w 
e 


\<»> 


(16) 


The  bracketed  term  in  (16)  corresponds  to  the  frequency  dependent  weight 
synthesized  by  the  delay  line  in  the  k—  channel.  We  define,  then. 


\(w)  =*  2) 

n=l 


+j  (n-lH0w 

A,  e 
k,n 


(17) 


to  be  this  frequency  dependent  delay  line  synthesized  weight.  Eq.  (17)  defines 
a  truncated  Fourier  series  expansion  for  each  weight  A^(w)  over  the  interval 
1  -  FBW/2 <  w  <  1  +  FBW/2.  Ideally,  A^(w)  can  be  made  to  approximate  (w) 
as  well  as  desired  by  choosing  4>q  appropriately  and  letting  be  arbitrarily 
large.  It  is  well  known  from  the  theory  of  such  series  expansions  that  in 
order  to  be  able  to  approximate  an  arbitrary  function  over  this  band,  <|>q 
<  2ir/FBW.  Note  that  if  we  choose  -  2tt/FBW,  then  A^w)  is  periodic  with 
period  FBW.  Consequently  any  Fourier  sum  approximating  AQk(w)  must  exhibit  the 
familiar  Gibb's  phenomenon  ripple  at  the  ends  of  the  band,  and  leads  to  oscil¬ 
lating  approximations  to  Aq^(w)  within  the  band.  Typically  the  optimum  weights 
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Aq^(w)  required  depend  only  on  the  array  dimensions  (or  array  environment,  if 
multipath  arising  from  boundaries  external  to  the  array  is  present)  and  are 
generally  slowly  varying  relative  to  2ir/FBW.  Consequently  it  is  normally 
advantageous  to  choose  <  2tt/FBW.  As  long  as  $q<  2tt/FBW,  it  is  clear  from 
the  theory  of  Fourier  series  expansions  that  the  approximation  of  A^(k)  to 
Aok(w)  will  improve  as  the  number  of  taps  increases.  It  is  the  intention  of 
the  following  sections  to  quantify  this  improvement  with  the  objective  of  using 
only  the  minimum  number  of  taps  for  a  specified  cancellation  performance. 

To  obtain  a  feel  for  how  well  the  delay-line  synthesized  weights  defined 
by  Eq.  (17)  approximate  the  ideal  weights,  A^(w),  consider  an  N-element  array 
and  a  single  interference  source  incident  from  an  angle  0,<f>.  Assume  = 

[1,  0,  .4. ,  0]  corresponding  to  turning  on  only  a  single  element  (i.e.,  earth 
coverage  for  an  array  on  a  satellite) .  Choose  the  array  phase  reference 
relative  to  y..;  then 
» 

A0i(w)  =  1 

i  j(\-Vn)w 

A0k(w)  “  ■  STT  e  »  k=2’  •••»  (18) 

By  comparing  Eqs.  (17)  and  (18),  several  comments  concerning  the  choice  of 
can  be  made: 

1)  Note  that  using  Nt  -  two  taps/element,  no  choice  of  (Jig  in  Eq.  (17) 

will  yield  exact  synthesis  of  (18)  unless  N-2.  In  this  case,  for  the  two 
element  array,  let  <)> q  =  Y2_YL  =  then  the  weighting  A^  2  *  0,  ^  =  0 

results  in  perfect  cancellation  over  FBW.  If  the  angle  of  incidence  defined 
by  9 , 4>  is  allowed  to  be  variable  over  the  FOV,  perfect  cancellation  of  the 
interference  waveform  will  be  achieved  only  for  selected  angles  of  incidence 
for  a  fixed  <£>q. 

2)  For  N-elements,  choosing  Nt=N  will  not  necessarily  yield  exact 
synthesis  of  the  optimum  weights  using  uniform  tap  spacing.  Conversely,  if 

a  different  tap-spacing  =  u8®d  ^or  each  element,  then  a  total  of 
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Nt=N  taps  are  sufficient  to  perfectly  cancel  the  interference  over  the  band 
for  selected  angles  of  incidence. 

3)  Since  any  angle  of  incidence  over  the  FOV  is  assumed  equally  probable, 
it  is  the  average  cancellation  which  can  be  achieved  over  this  range  of  angles 
which  becomes  important.  In  this  case,  the  use  of  uniform  tap-spacing  does 
not  appear  to  be  overly  restrictive. 

4)  Referring  to  (17)  and  (18),  an  interesting  graphical  aid  in  visualizing 
how  the  choice  of  affects  the  synthesis  of  the  optimum  weights  can  be  devel¬ 
oped.  Consider  a  polar  plot  (amplitude/phase)  of  Eqs.  (17)  and  (18),  where 

for  simplicity  we  consider  only  the  case  of  2  and  3  taps/element.  Fig.  3a 
illustrates  the  graphical  synthesis  of  the  A^^(w)  using  2  taps/element,  and 
Fig.  3b  using  3  taps/element.  The  optimum  weight  A^(w)  has  a  total  phase 
traversal  of  (-^“Y^)  *  FBW  radians.  For  A^  ^  and  A^  2  4  0,  as  would  generally 
be  the  case  for  a  fixed  then,  at  best,  the  vector  sum  of  the  delay  line 
phasors  can  only  approximate  AQk(w).  The  two-element,  two-tap  special  case, 
with  4>q  -  Y2“Y1»  becomes  clear  from  Fig.  3a,  in  which  case  choosing  A^  ^  =  0 
results  in  an  exact  synthesis  of  Aq^w) •  Fig-  3b  for  the  three  tap  case 
illustrates  the  added  degree  of  freedom  obtained  by  adding  the  additional  delay 
line.  In  both  examples,  has  been  chosen  too  large,  and  the  figures  show  the 
resultant  too- rapid  phase  variation  which  would  be  obtained.  This  phase- 
amplitude  plot  will  be  used  in  the  later  sections  to  illustrate  the  performance 
of  the  delay  line  processor  relative  to  optimum  weighting  for  specific  array 
antennas. 
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optimum  weight  Aq(w) 
/Phase  plot  using  tw 


Fig.  3.  Synthesis  of  the 
delay  lines:  (a)  Amplitude 
(b)  Three  taps/element. 


III.  THE  TWO-ELEMENT  ARRAY 


The  two-element  array  presents  the  simplest  array  antenna  configuration 
which  can  be  used  to  achieve  radiation  pattern  nulling.  Though  simple  in 
concept,  a  detailed  analysis  of  its  performance  characteristics  provides  valu¬ 
able  insight  into  the  mechanism  of  delay  line  frequency  compensation.  Prac¬ 
tical  considerations  dictate  that  only  the  minimum  number  of  taps  required  to 
achieve  the  desired  array  performance  be  used.  With  this  goal  in  mind,  the 
two-element,  two-tap  problem  will  be  studied  in  some  detail.  For  this  simple 
case,  analytical  solutions  for  the  delay  line  synthesized  weights  can  be 
obtained  together  with  analytical  solutions  for  the  optimum  weights.  In 
comparing  the  two  solutions,  it  will  become  clear  that  exact  synthesis  of  the 
optimum  weights  via  two  taps  on  a  delay  line  is  possible  when  the  tap  spacing 
is  matched  to  the  wavefront  angle  of  arrival.  However,  for  the  practical 
case,  this  angle  of  arrival  is  not  a  priori  known,  so  that  this  matched  tap 
spacing  could  not  be  used.  In  practice,  a  fixed  delay  line  spacing  must  be 
implemented  to  compensate  for  an  arbitrary  angle  of  arrival.  For  our  purposes, 
we  define  an  "optimum"  delay-line  tap  spacing  relative  to  this  latter  crite¬ 
rion.  That  is,  the  optimum  delay-line  tap  spacing  minimizes  the  average 
cancellation  as  a  function  of  signal  wavefronts  Incident  from  anywhere  over 
the  antenna  FOV.  Solutions  for  the  optimum  tap-spacing  satisfying  this  crite¬ 
rion  will  be  obtained  for  the  two-element  array  for  two  taps  on  the  delay  line, 
and  extrapolated  to  the  N-element  array  employing  a  delay  line  with  Nt  taps 
via  physical  considerations.  It  is  also  possible  to  obtain  analytical  solu¬ 
tions  for  the  eigenvectors  and  eigenvalues  of  the  correlation  matrix  defined 
at  the  tap  outputs  for  the  two-element,  two-tap  array.  Knowledge  of  these 
eigenvectors  and  eigenvalues  allows  a  detailed  examination  of  the  adapted 
weights  to  be  made  on  the  basis  of  Eq.  (13).  In  particular,  the  dependence  of 
the  array  cancellation  performance  on  processor  dynamic  range  discussed  in  the 
previous  sections  can  be  developed.  Finally,  by  extending  the  two-tap  results 
to  the  more  general  N  >  2,  a  quantitative  assessment  of  the  improvement  in 
cancellation  realized  by  Increasing  the  number  of  taps  can  be  made. 
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In  the  following,  we  first  develop  the  characteristics  of  the  eigenvalues 
and  eigenvectors  of  R  characterizing  the  performance  of  the  two-element,  two- 
tap  topology,  and  examine  each  respective  eigenvector  in  some  detail.  As 
might  be  expected,  the  eigenvector  e^,  corresponding  to  the  lowest  eigenvalue, 
contains  the  basic  information  on  how  to  form  nulls  with  delay  line  compen¬ 
sation.  The  transition  of  the  adapted  weight  vector  A  defined  in  Eq.  (13) 
to  the  optimum  weighting  defined  by  e^  is  then  considered  as  a  function  of 
processor  dynamic  range,  ^ We  show  that  as  °°»  A  -*■  e^  to  within  a 

constant  factor.  Finally,  given  the  case  us*,.,,  >>  1,  we  determine  the  choice 

MAX 

of  <j>Q  which  maximizes  the  average  interference  cancellation  for  sources  inci¬ 
dent  from  anywhere  over  a  fixed  FOV. 

A.  Eigenvectors  and  Eigenvalues  of  R  for  a  Two-Tap  Delay  Line 


Consider  the  two-element,  two  tap  topology  illustrated  in  Fig.  4.  For 
convenience  we  define  Eg  =  v/G^R^  and  Ay  =  Y^  -  Y2*  The  solution  is  simplified 
considerably  with  little  loss  in  generality  if  we  ignore  the  effects  of  thermal 
noise.  For  a  single  interference  source  incident  on  the  array,  the  4x4 
correlation  matrix  characterizing  the  delay  line  outputs  is  given  by  (see  Eqs. 
11  and  12) 


=4  "  (E0 


)/2  <  e 


+j  (  Yk-  Y )w  -j  (n-m)<J>0w 


(19) 


where  we  have  defined  v  ,  =  Y.  and  y  ,  =  Y  •  If  the  eigenvectors  and  eigen- 

K  q,i  q 

values  of  R^  can  be  determined,  then  using  Eq.  (13)  the  complete  solution 
for  the  adapted  weights,  and  other  physically  observable  parameters,  can  be 
obtained.  An  exact  analytical  solution  for  these  four  eigenvalues  and  eigen¬ 
vectors  is  difficult.  However,  they  can  be  estimated  quite  accurately  using  a 
perturbation  expansion  of  R^  on  FBW. 

To  accomplish  this,  we  expand  the  integrand  of  (19)  about  w*l.  Define 
Rq  according  to 
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Fig.  4.  Delay  line  topology  and  notation  for  the  two-element 
two  taps/element,  delay  compensated  array. 


R~  s  R 
=0  =4 


FBW=0 


(20) 


Then  R^  can  be  written  In  the  series  expansion 


R.  =  R.  +  AR.  +  AR,  +  AR,  =  R~  +  AR 

=4  =0  =2  =4  =6  =0  = 


(21) 


2  4 

where  AR  is  proportional  to  FBW  ,  AR  to  FBW  ,  and  AR, 
— L  —4  — D 

powers  of  FBW  are  not  present  since  <Aw  >  =  0  when  n  is 
The  eigenvalues  and  eigenvectors  of  R^  can  be^ estimated 


to  FBW**.  (Note,  odd 
odd ,  where  Aw=w-1 . ) 
by  the  expressions^ 


s 


i 


(22a) 


i  =  1,2, 3, 4 


^i  = 


e,  n  +  2  - — \ -  (e.  +  „  *  AE  *  e,  n)  e 

-i,0  s  ~s  n  “J»°  -i*°  ~J»° 

j=l  1 


(22b) 


where  the  s^  q  and  e^  i=l,  ...,  4  are  the  solutions  to  the  unperturbed 
problem,  (i.e.,  the  eigenvectors  and  eigenvalues  of  Rq)  and  the  and  e^, 
i“l,  . ..,  4  are  the  estimated  eigenvalues  and  eigenvectors  of  R.  Note  by  the 
form  of  AR  that  these  estimates  are  correct  to  the  order  of  FBW  .  The  analysis 
leading  up  to  the  expression  for  the  s^  and  e^  is  long  and  tedious.  The 
procedure  is  outlined  in  Appendix  B.  The  eigenvalues  take  the  form 


s2s  EQ2  (FBW  •  Ay)2  (1  +  <$>Q2/Ay2)/24 


(23) 
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23 


(26) 


S10  "  S1  +  s2  +  S3  +  S4 

As  FBW  increases,  since  s1Q  is  independent  of  FBW,  generally  s^  begins  to 
decrease  and  82*  a ^  and  s^  Increase  so  that  (26)  is  satisfied.  Hence  since 
s^  is  the  smallest,  the  approximation  to  s^  in  (23)  is  the  most  accurate. 

We  now  show  that,  as  contrasted  to  the  eigenvalues  of  gA  (where  gA  is  the 
N  x  N  correlation  matrix  defined  at  the  antenna  output  ports),  the  eigenvalues 
of  R^  cannot  be  associated  directly  with  output  power  relative  to  thermal 
noise.  To  see  this,  we  note  by  definition  that 


whereas  the  output  power  relative  to  thermal  noise,  (i/N)^,  when  e^  is  applied 
as  a  weight  is  given  by 

e  +*  R,  •  e, 

(i/N)  -  - - — - -  ,  (28) 

<A  (w)  *  A(w)> 


where  A(w)  is  the  delay  line  synthesized  weight  defined  in  Eq.  (17).  Thus 
(I/N)^  is  given  by,  using  Eq.  (8), 


(i/N) ±  = 


Sa 


(29) 


5n  ’  % 


and  R^  is  not  a  diagonal  matrix.  Consequently 


(X/N)±  =  Si 


t 

%  *  ^i 


*  *N  *  «i 


(30) 
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Hence  although  (I/N)^  Is  proportional  to  s^,  a  strict  equality  does  not  hold. 
Consequently,  we  Interpret  the  s^  relative  to  output  Interference/thermal  noise 
ratio  in  the  loose  sense  indicated  by  Eq.  (30). 

B.  Interpretation  of  the  Eigenvectors  and  Eigenvalues 

We  have  seen  from  Eq.  (13)  that  the  adapted  weight  vector  which  minimizes 
the  array  output  power  is  comprised  of  a  sum  of  eigenvectors,  suitably  weighted. 
It  is  instructive  to  examine  the  delay  line  synthesized  weight  and  resultant 
output  power  when  each  eigenvector,  e^,  i=l»  •••»  4,  in  turn,  is  used  as  a 
weight  vector.  The  results  give  some  insight  into  the  physical  significance 
of  the  eigenvectors. 

1.  Tapped  Output  Weighted  by  e^  and  e^ 

The  first  two  eigenvalues,  s^  and  s^,  can  be  interpreted  similarly  to  the 

first  two  eigenvalues  of  the  correlation  matrix  defined  at  the  antenna 

element  outputs.  That  is,  we  note  that  s^  2  E^  represents  the  maximum 

output  power  obtainable  from  the  array.  This  output  results  from  effective 

weights  (w)  and  A2(w)  which  cause  the  two  weighted  channels  to  add  in  phase 

at  the  output.  This  results  in  a  radiation  pattern  having  a  beam  maximum 

pointed  toward  the  interference,  which  is  relatively  independent  of  frequency. 

2 

s_,  represents  the  familiar  array  dispersion  phenomenon,  which  varies  as  FBW 
compounded  by  dispersion  Introduced  by  the  delay  lines.  It  is  zero  for  zero 
bandwidth  and  increases  as  the  square  of  the  bandwidth.  Applying  e^  as  a 
weight  vector  in  essence  puts  a  "zero-bandwidth"  null  on  the  source  and 
generates  a  "monopulse"  type  radiation  pattern.  Examination  of  A1 (w)  and 
A2(w)  using  e^  as  a  weight  at  the  tap  outputs  shows  that  A^(w)  and  A2(w)  are 

only  weakly  dependent  on  w.  Consequently,  weighting  characterized  by  e_^  and 

t  J 

e^  represent  the  most  interesting  cases  relative  to  broadband  nulling. 

Furthermore,  since  s^,  s2  »s^  s^,  it  is  only  e ^  and  e^  which  contributed 

significantly  to  the  adapted  weights.  (Unless,  of  course,  -*■  0,  in  which 

t  +  U 

case  _e^  •  V  and  •  V  +  0  in  (13);  the  solution  for  A  then  reduces  to  that  of 

frequency  independent  weighting  at  the  array  output  for  which  s2  characterizes 

the  adapted  array  output  power.  This  case  is  well  understood.) 


25 


2.  Tapped  Output  Weighted  by  e^ 

Using  the  expression  for  e^  from  Eq.  (24)  in  Eq.  (17)  results  in  the 
delay  line  synthesized  weights  A^(w)  and  A^(w)  given  by: 

.  *,  .  1  -jA  Yf  -^0(wtal) 

A1  (w)  *  2  e  (e  "  ^ 


A2*(w)  “  j  [1  -  e 


-J4>0  (w-1) 


Ola) 

(31b) 


Equation  (31)  leads  to  an  output  voltage  vs.  frequency  relationship,  V^(w), 
given  by 


V(w)  =  A1  (w)  E^(w)  +  A2  (w)  E2(v) 

.  jy2  j  Y.(w-l)  jY2(w-l) 

»  i  e  [e  ~  1]  [e  -  e  ] 


(32) 


+ 

Hence  when  e3  is  used  as  weight  vector,  the  resultant  output  at  w-1  is  zero 

by  virtue  of  the  fact  that  each  Ak(w»l)*0,  i.e.,  the  array  turns  off  at  band 

center,  as  opposed  to  cancelling  the  weighted  interference  between  channels. 

Thus  it  would  be  undesirable  to  use  e^  to  synthesize  the  weights  since  the 

user  signals  as  well  as  the  Interference  signals  are  nulled.  In  order  to 

assure  that  the  adapted  solution,  Eq.  (17)  employs  e_  in  the  solution  as 

+ 

little  as  possible,  we  choose  V  so  that  V  •  e3»0  which  from  Eq.  (13)  prevents 
e^  from  contributing  to  the  adapted  weight.  This  condition  will  be  satisfied 
if  we  choose  V  of  the  form: 


V+ct 


(33) 


where  and  v2  are  arbitrary  complex  constants.  Physically  this  choice  of 
V  corresponds  to  a  quiescent  weighting  behind  a  specific  element  so  that  the 
ouputs  of  each  tap  add  in  phase  in  the  absence  of  interference. 


For  example,  consider  a  wavefront  normally  incident  on  the  array.  Two 
possible  choices  of  V+  can  be  used  to  turn  on  array  element  number  one,  if  an 
earth  coverage  quiescent  radiation  pattern  is  desired.  One  choice  would  be: 

V  =  [1,0, 0,0],  which  is  inconsistent  with  Eq.  (33)  above.  For  this  choice  of 

V*,  an  allowable  adapted  weight  vector  A  which  results  in  minimum  output  power 

t  j4>n 

would  be  A  *  [1,0, -e  u f  0],  This  adapted  weight  turns  off  the  array  at  band 

•f*  *J* 

center  (i-e..  A, (w*l)=0  and  Ao(w=l)=0).  However,  choosing  V  of  the  form  V  = 
[1,0,  +  e  °,0]  prevents  this  situation  from  occurring,  since  A  can  no  longer 
contain  a  component  of  e^*  The  latter  choice  of  assures  V+*  e^  =  0. 
Additionally,  this  choice  of  V  gives  maximum  output  of  the  delay  lines  in  the 
quiescent  mode  of  operation. 

When  e~  is  applied  as  a  weight,  s,  governs  the  adapted  output  power  of 
J  ?  2  4 

the  array.  Note  the  dependence  s^  «  4>q  A -y  FBW  .  As  c|)q  "*•  0,  then  s^  -*■  0. 

Thus,  for  *  0,  the  only  possible  means  of  obtaining  zero  output  power  when 
FBW  f  0  is  for  the  array  to  turn  itself  off  by  adjusting  the  output  of  the 
taps  (the  delay  line  is  only  a  simple  power  divider  in  this  limiting  case)  to 
be  180°  out  of  phase  with  each  other. 

3.  Tapped  Output  Weighted  by  e^ 

Examination  of  the  eigenvalue  s^  shows  that  it  goes  to  zero  when  4>q  =  0 
and  <j>Q  ■  Ay,  Independent  of  FBW.  The  case  s^  =  0  when  4>q  *  0  is  similar  to 
the  situation  discussed  above.  It  leads  to  the  trivial  solution  of  zero 
output  power  results  from  turning  off  the  array.  The  case  s^  =  0  when  * 

Ay  leads  to  the  desired  optimum  result.  That  is,  applying  e^  as  weight  leads 
to  the  optimum  A^(w)  and  synthesized  using  two  taps/element.  These  are 

given  by 

*  e-j4y  -j*0(v-l) 

Ax  (w)  *  . . . . t  (1  -  4>qMy)  -  (1  +  <t>0/&Y)e  ] 

2\/l  +  t02/4>2 

(34) 


A-  (w) 


2\4  + 


[(1  +  i >n/Ay)  -  (1  - 


2,.  2 

i  My 


<b0/A  y)e 


-J4>0(w-1) 
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Comparing  A^w)  and  A2(w)  in  Eq.  (34)  to  the  optimum  A^Cw)  and  Aq2(w)  which 
null  the  interference  sources  over  the  entire  band  yields  exact  comparison 
when  ■  | Ay | -  Note  that  since  Ay  “  kgD  sin  then  the  delay  line  synthe¬ 
sized  weighting  is  independent  of  the  sign  of  0.  Hence  when  Ay  >  0  and  4>q  * 
Ay,  Eq.  (34)  reduces  to 


A1  (w) 


1  e-jAyw 


A2  (w) 


and  when  Ay  <  0  and  4>  q 


|  Ay |,  we  obtain 


Aj* (w) 


=  -i-  e 
v/r 


e-JAY 


4-jAyw 

6 


(35) 


(36) 


Substituting  either  (35)  or  (36)  into  Eq.  (16)  yields  perfect  cancellation  of 
the  interference  sources  over  the  entire  band. 

4.  Simulations 

The  approximations  to  the  Sj^  and  e^  defined  in  Eqs.  (23)  and  (24)  have 

been  verified  by  computer  simulation  using  the  exact  solution  outlined  In 

Section  II  for  the  N-element-Nt  tap  topology.  Figure  5  Illustrates  the  general 

variation  of  the  eigenvalues  as  4>q  Is  varied  for  parameter  AyFBW  =  1.9. 

Increasing  AyFBW  leads  to  a  much  more  pronounced  resonance  for  the  lowest 

eigenvalue  s, .  Note  that  when  <f>_  ■  |  Ay  |  ,  the  resultant  output  power  *  0, 

4  t  u 

as  8^  *  0  if  e^  Is  applied  as  weights  to  the  tapped  delay  line  outputs.  For 
<1>q  >>  I  Ay  |,  s^  a  FBW^  as  predicted  by  Eq.  (23). 
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Fig.  5.  Eigenvalue  characterization  vs  tap  spacing  for  the 
two-element  array  using  two  taps /element.  Ly  •  FBW  =  1.9. 
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C.  Algorithm  Dependence/Performance  Evaluation 


According  to  the  above  discussion,  for  a  specified  positioning  of  the 

interference  source  it  is  clear  that  optimum  weighting  results  if  =  |Ay| 

t  U 

and  the  eigenvector  e^  is  applied  as  a  weight  to  the  tapped  delay  line  outputs. 

However,  in  practice,  the  position  of  the  interference  source  is  not  apriori 

known.  Furthermore,  it  is  desired  that  the  weights  be  set  autonomously 

according  to  some  prescribed  algorithm.  Both  of  these  factors  dictate  that 

perfect  cancellation  as  governed  by  e^  for  q  =  I  Ay |  will  not  generally  be 

achieved.  For  our  purposes,  we  have  chosen  to  fix  the  delay  line  phase  <t>g 

according  to  some  prescribed  value  (to  be  determined  later) ,  and  set  the  weights 

according  to  the  Applebaum-Howells  criterion  specified  by  Eq.  (6).  Consequently 

the  adapted  weights  are  proportional  to  the  sum  of  all  the  eigenvectors,  each 

weighted  according  to  l/Q+ps^).  Generally  s1>>s2»  s^  and  s^,  and  we  choose 

V  according  to  Eq.  (33),  so  that  e^+*  V  =  0.  Then  the  adapted  weights  take 

the  form 


A  = 


1  +  ps„ 


(£2' 


V)  e,  + 


1  +  MS. 


<2* 


D 


(37) 


The  solution  governed  by  Eq.  (37)  is  broadband  only  when  e,  dominates.  For 

f 

this  to  be  true  we  require  ms^  >>  1  and  <f>^  large  enough  so  that  •  V  is  not 
small  (note,  from  Eqs.  (24)  and  (34),  |e^  •  vj  «  (<j>g/Ay)  ,  which  approaches 
zero  as  4>q  -*■  0).  In  turn,  the  degree  to  which  ps2  is  much  greater  than  unity 
depends  on  the  loop  gain  m  (the  value  of  which  is  more  conveniently  specified 
in  terms  of  the  processor  dynamic  range,  Ms^^),  and  on  the  bandwidth  FBW. 
Consequently  the  performance  of  the  delay  line  processor  used  in  conjunction 
with  an  Applebaum-Howells  power  inversion  algorithm  can  be  expected  to  be  a 
function  of  the  variables  Ms^^  ,  4>qMy  and  FBW.  Since  the  eigenvalues  them¬ 
selves  are  a  function  of  FBW  by  way  of  the  parameter  Ay* FBW,  this  latter 
parameter  offers  a  more  universal  characterization  of  the  results. 

The  above  discussion  illustrates  that  the  cancellation  C  is  a  function 
of  the  three  parameters  MSmax,  4>q/Ay  and  Ay* FBW,  where  C,  defined  as 
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c 


(38) 


a/N)a 
=  <I/N)b 


is  the  ratio  of  the  interference  to 
interference  sources.  In  Eq.  (38), 


«/N)b 


thermal  noise  before  and  after  nulling  the 

(I/N),  and  (I/N)  are  given  by 
d  a 

A+*  R.*  A 

(I/N)  =  ~r~— -  (39) 


It  is  possible  to  then  evaluate  Eq.  (38)  analytically  using  Eq.  (37)  for  A. 
However,  the  result  is  tedious  and  will  not  be  carried  out  here.  Rather  we 
illustrate  in  Figs.  6,  7  and  8  the  dependence  of  C  on  these  parameters  as  each 
is  varied  for  fixed  values  of  the  remaining  ones.  The  results  were  obtained 
from  computer  simulation,  and  are  not  dependent  on  the  perturbation  expansion 
technique . 

Figs.  6a  and  b  illustrate  the  dependence  of  C  as  $Q/Ay  is  varied  for 
Ay  •  FBW  =*  1.0  and  4.0,  respectively,  with  as  a  parameter.  Note  that,  as 

predicted  from  the  behavior  of  e^,  when  Ay  •  FBW  is  large,  and  >>  1,  a 

pronounced  resonance  in  the  dependence  of  C  on  ^/Ay,  indicating  broadband 
cancellation  is  achieved.  The  width  of  the  resonance  is  dependent  on  the  values 
of  Ay  •  FBW,  as  can  be  seen  by  comparing  Figs.  6a  and  b.  Note  that  for  Ay  •  FBW 
<  1,  the  resonance  is  quite  broad,  indicating  that  employing  two  taps/element 
is  adequate  to  give  good  cancellation  over  a  wide  range  of  interference- 
source- locations,  with  fixed  tapped  spacing.  However,  as  the  aperture  disper¬ 
sion  increases  (either  by  increasing  FBW  or  —  sin  0),  more  than  two  taps/ 
element  would  be  required  to  achieve  deep,  broadband  cancellation  over  a  wide 
FOV. 

Figs.  7  and  8  illustrate  the  same  general  behavior  as  Fig.  6,  only 
presented  using  differing  values  of  the  three  parameters.  Fig.  7  illustrates 
the  dependence  of  C  vs  $q/Ay  as  the  parameter  Ay  •  FBW  changes,  for  a  fixed 
power  level  whereas  Fig.  8  shows  the  dependence  of  C  on  with 
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Fig.  6.  Cancellation  (dB)  vs  tap  spacing  for  the  two-element 
two  tap  array  with  ps.,.,,  and  Ay  ♦  FBW  fixed:  (a)  Ay  •  FBW  *  1 


Fig.  7.  Cancellation  (dB)  vs  tap  spacing  for  the  two-element 
two  tap  array  with  ysWJV  and  Ay  *  FBW  fixed:  (a)  psUAV  ■  25  d 


S  10  IS  to  *•  to  M  40  *» 

mW«> - - 


Fig.  8.  Cancellation  (dB)  vs  ^or  C^e  two-element«  two  taP 

array  with  Ay  •  FBW  and  tap  spacing  fixed:  (a)  Ay  •  FBW  ■  1.0, 
(b)  Ay  •  FBW  =  4.0. 
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4>  /A  y  ^9  a  parameter  and  Ay  •  FBW  fixed.  Note  from  Fig.  8  that  a  limit  exists 
where  C  becomes  independent  of  power  level.  This  limit  is  a  fundamental  limit 
based  on  the  array  and  delay-line  geometry,  on  the  bandwidth,  and  on  the 
interference  source  location.  It  could  not  be  improved  by  use  of  another 
algorithm.  In  this  region,  one  is  in  essence  applying  the  eigenvector  e^*  as 
weights  to  the  tap  outputs.  Fig.  8  illustrates  that  if  the  processor  gain  y 
is  set  so  that,  for  example,  >30-35  dB,  then  AO  dB  cancellation  is 

possible  on  these  larger  sources  over  a  broad  range  of  <I>q/A  y,  corresponding 
physically  to  a  wide  range  of  incidence  angles.  This  indicates  that,  in 
general,  a  set  of  fixed  delay  line  parameters  exist  where  significant  cancel¬ 
lation  can  be  achieved  compatible  with  a  wide  range  of  source  scenarios. 
Optimization  of  this  choice  of  delay  line  parameters  for  differing  source 
scenarios  is  the  objective  of  the  following  section. 

D.  Optimize  <j>Q  for  a  Predetermined  FOV 

The  above  results  have  indicated  that  if  A y  were  known  (i.e.,  the  inci¬ 
dent  angle  of  arrival),  then  choosing  ♦q  “  | A  yI  yields  optimum  wideband 
cancellation  for  the  two  element  array,  using  two  taps/element  (Nt  >  2  will  be 
considered  in  the  next  section).  Generally,  however,  the  interference  may  be 
located  anywhere  over  the  antenna  FOV  and  the  angle  of  arrival  is  not  a  priori 
known.  Hence  it  is  of  interest  to  optimize  the  choice  of  4>q  for  best  average 
cancellation  over  a  range  of  Ay,  say  -A  Ym<  Ay^Ay^.  Since  C  is  a  function 
of  VJSMAX  aa  well  as  A  y  *  FBW  and  4>qMy*  optimization  on  is  generally  compli¬ 
cated.  However,  in  the  limit  >>1,  then  the  adapted  weight  is  approxi¬ 

mately  the  eigenvector  e^.  In  this  case,  substituting  the  expression  for  e^ 
into  Eq.  (28),  the  interference/thermal  noise  ratio  after  adaption  can  be 
shown  to  be 

1  +  *02/Ay2 

<I/N)a  =  S4  ■  ~ 2°  2  *  (40> 

4>0  My 
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where  we  have  approximated  <  (w)*A(w)  >  as  A** A  evaluated  at  w=l,  and  is 
defined  in  Eq.  (23).  It  is  straightforward  to  average  (I/N)  over  -Ay  <  Ay 

3l  m 

<  Aym  and  maximize  this  average  with  respect  to  $q.  This  result  is  given  by 

=  0. 775Ay  .  (41) 

u  m 

Hence  best  average  cancellation  assuming  the  interference  sources  located  any¬ 
where  in  the  range  -0  <  0  <0  results  from  choosing  —  three  quarters  of 

mm  u 

Ay^.  Note  since  C  is  independent  of  the  sign  of  Ay,  only  a  single  delay  line 
is  required  to  provide  cancellation  over  both  positive  and  negatives  values 
of  0. 

Fig.  9  illustrates  the  resultant  cancellation  vs  0  (expressed  in  terms  of 

Ay/Ay  )  using  one  and  two  taps/element  as  the  incidence  angle  is  varied  over 
m 

the  interval  -0  <  0  <  0  .  AyFBW  is  fixed  at  6.0  and  <)>_  is  fixed  at  0.775 
mm!  0 

Ay  .  Only  positive  values  of  0  are  plotted  since  the  results  are  symmetric 
m 

about  0=0.  The  choice  of  <J>Q  used  for  two  taps/element  results  in  the  best 
average  cancellation  over  the  FOV.  Worst  case  results  occur  when  0=0^. 

Since  the  optimization  criterion  leading  to  4>q=0.775  Ay^  was  somewhat  arbi¬ 
trary,  one  might  consider  increasing  <t>Q  and  minimizing  the  peak  cancellation 

over  -0  <  0  <  0  .  The  resultant  <t>«  would  be  greater  than  0.775  0  . 
m  m  U  m 

E.  Extensions  to  >  2 

Depending  on  the  parameter  Ay  *FBW  and  the  cancellation  required  for  a 

m 

particular  system,  using  only  two  taps/element  might  not  yield  adequate  cancel¬ 
lation.  In  this  section  we  briefly  consider  using  multiple  taps/element  for 
the  two-element  array.  It  will  become  clear  later  that  the  general  trends  and 
baseline  parameters  characterizing  the  two-element  array  are  directly  appli¬ 
cable  to  the  N-element  array. 

Consider,  then,  the  performance  of  the  two  element  array  when  A^(w)  and 
A2<w)  are  synthesized  using  taps/element,  and  assume  the  tap  spacing  4>q 
between  taps  is  uniform.  It  is  clear  from  Fig.  (3)  that  increasing  Nt  leads 
to  a  better  approximation  of  Aq^(w)  by  A^(w).  Of  course,  if  4>0  =  I  I , 
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Fig.  9.  Cancellation  (dB)  vs  angle  of  arrival  for  the  two- 
element  array  using  frequency  independent  weighting  and  delay 
line  weighting  with  two  and  three  taps  per  element. 
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there  is  no  need  for  Nf  <  2  as  the  optlmim  weights  are  synthesized  exactly  for 

this  case.  For  this  reason  we  fix  the  delay  line  spacing  <t>n  <  Ay  and  vary 

u  tn 

0  over  the  FOV  as  defined  above.  For  N  =3,  we  choose  *  tt(0.775  Ay  )  = 

t  u  z  m 

0.3875  Ay^  —  i.e.,  one  half  the  optimum  value  of  for  Nt=2.  This  approxi¬ 
mates  a  nearly  uniform  separation  in  when  viewed  as  A  y  varies  over  the 
field  of  view.  This  choice  of  <1>q  now  yields  perfect  cancellation  for  two 

positive  incidence  angles:  Q  =  0.775  Q  and  0  *  0.3875  0  .  The  cancellation 

m  m 

realized  using  N^=3  as  0  varies  over  the  field  of  view  is  illustrated  in  Fig, 

9  for  comparison  with  the  case  Nt=2  discussed  above.  Note  the  marked  improve¬ 
ment  obtained.  This  improvement  results  from  the  additional  freedom  to  cancel 
the  interference  source  exactly  at  a  second  positive  angle  of  arrival  over  the 
FOV,  which  reduces  the  average  cancellation  over  the  FOV  considerably. 

Because  of  our  choice  of  4>q,  worst  case  cancellation  again  occurs  at  0=0^. 

Of  course,  by  adjusting  <|>q  appropriately  the  dispersion  at  0=0^  can  be  reduced, 
at  the  cost  of  degraded  performance  towards  the  center  of  the  FOV. 

Noting  that  worst  case  cancellation  for  our  choices  of  <t>g  occur  when 

0=0^,  the  cancellation  obtained  for  this  angle  of  arrival  can  be  used  to 

characterize  quantitatively  the  cancellation  dependence  on  the  number  of  taps 

as  Ay  • FBW  is  varied.  To  this  end.  Fig.  10  illustrates  cancellation  as 
m 

A^*FBW  is  varied  using  N  *1,  2  and  3  and  the  same  choices  for  as  in 

Fig.  9  (note,  the  curves  for  3  and  6  elements  will  be  discussed  later).  To 

obtain  the  result  illustrated  for  the  two  element  array,  the  angle  of  arrival 

0=0  and  aperture  diameter  D/X  is  fixed,  and  FBW  is  varied  accordingly.  The 
m  D 

utility  of  the  parameter  2ir •—  sin  0  *FBW  in  characterizing  the  delay  line 

a  in 

performance  becomes  clear  from  the  results  shown  in  the  figure.  For  example, 

assume  the  antenna/delay  line  processor  is  to  be  designed  to  achieve  30  dB 

cancellation  for  interference  sources  located  anywhere  over  the  FOV.  The 

system  requirements  (i.e.,  signal-interference  resolution  required,  system 

bandwidth,  etc.)  determine  the  value  of  Ay  *FBW  required  to  realize  a  given 

m 

performance  specification.  Clearly  if  A y  • FBW  <  0. 25,  then  frequency  indep¬ 
endent  weights  would  suffice  to  realize  the  desired  cancellation.  If  0.25 
<  Ay^'FBW  <  4,  then  a  delay  line  using  two  taps/element  should  be  employed; 
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if  4  <  A  y^'  FBW  <9,  three  taps/element  should  be  used,  and  so  forth.  The 
results  for  the  three  and  six  element  arrays  shown  in  Fig.  10  (to  be  consid¬ 
ered  later)  Indicate  this  characterization  to  be  somewhat  invariant  to  the 
number  of  elements. 

One  final  observation  can  be  made  concerning  the  delay  line  performance 

results  for  N  >  2;  namely,  the  cancellation  varies  as  a  function  of  (A y  • 

£  t 

FBW)  ,  where  t  is  approximately  given  by  2(2N  -1).  Thus  for  frequency  indep- 

2  ^ 

endent  weights  C  «  (Ay  •  FBW)  ;  for  delay  lines  using  two  taps/element,  C  « 

6  ® 

(Ay^'FBW)  ,  and  so  forth  for  Nt  >  2.  This  dependence  can  also  be  developed 

analytically  using  a  bandwidth  expansion  for  the  correlation  matrix  similar  to 

that  used  in  Eq.  (26).  Generally,  for  N-l  sources,  R  has  N*Nt  eigenvalues. 

The  nulled  output  of  the  delay  line  processor  is  generally  characterized  by 

the  smallest  eigenvalue,  s  ,  which  it  can  be  shown,  for  the  two-element 

2  (2N 

array  varies  according  to  (A^'FBW)  t  ,  leading  to  the  above-mentioned 
cancellation  dependence.  Knowledge  of  this  variation  with  Ay  •  FBW  is  often 
quite  useful  in  scaling  a  known  performance  from  one  aperture  size  or  band¬ 
width  to  another  as  a  function  of  the  number  of  taps. 


IV.  THE  N-ELEMENT  ARRAY 


The  results  of  the  previous  section  yield  quantitative  estimates  of  the 
performance  characteristics  of  a  two-element  array  using  either  a  1,  2  or  3 
tapped  delay  line  for  each  element.  When  the  number  of  array  elements  exceeds 
two,  it  is  difficult  to  obtain  analytical  results,  so  that  the  majority  of  the 
analysis  for  N  >  2  is  based  on  numerical  simulations.  When  N  exceeds  two,  many 
configurations  for  array  element  positioning  become  possible,  and  it  is  intrac¬ 
table  to  consider  all  possibilities.  For  this  reason,  we  concentrate  our 
attention  on  the  three  classes  of  planar  arrays  illustrated  in  Fig.  11,  each 
consisting  of  only  a  relatively  small  number  of  elements.  For  reasons  which 
will  become  clear  later,  we  classify  the  arrays  according  to  their  "angular 
dispersiveness",  i.e.,  according  to  how  spread  out  the  N  elements  are  relative 
to  their  angular  distribution  over  a  two  dimensional  planar  aperture.  For 
example,  the  simplest  configuration  having  the  least  angular  spread  for  N- 
elements  distributed  over  a  planar  surface  would  be  the  equi-spaced  linear 
array  of  elements  illustrated  in  Fig.  11a.  One  might  view  such  an  array  as 
consisting  of  two  "radial  arms"  emanating  from  a  common  origin,  on  which  the 
elements  are  positioned.  The  second  configuration,  having  the  next  highest 
degree  of  angular  spread  has  three  radial  arms,  equi-spaced  120°  in  angular 
separation,  on  which  the  array  elements  are  positioned.  This  class  of  arrays, 
which  we  refer  to  as  triangular  arrays,  is  illustrated  in  Fig.  lib.  Finally, 
Fig.  11c  illustrates  four  other  array  types,  each  characterized  by  greater  than 
3  radial  arms  emanating  from  a  common  origin,  where  the  elements  are  positioned 
at  the  end-points  of  these  arms.  We  classify  these  as  "more  complex  arrays" 
relative  to  the  triangular  and  linear  arrays.  Carrying  through  this  delinea¬ 
tion  to  N  elements,  N  >>  1,  spread  out  uniformly  over  an  aperture  of  fixed 
diameter,  would  lead  to  N  radial  arms,  of  varying  radii,  characterizing  the 
planar  array. 

The  basic  idea  behind  this  categorization  of  the  array  types  considered  is 
as  follows:  Eq.  (15)  defines  the  optimum  weighting  for  an  arbitrary  N-element 
array.  When  the  array  elements  are  positioned  according  to  simple  geometrical 
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Fig.  11.  Three  classes  of  planar  arrays  for  use  with  delay  line 
compensation  evaluation:  (a)  Equi-Spaced  Linear  Arrays,  (b)  Triangular 
Arrays,  (c)  More  complex  arrays. 
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relationships  one  would  anticipate  that  the  frequency  variation  of  the  optimum 
weights  would  be  simply  related  between  element  outputs.  (See,  for  example, 

Eq.  (18)  specifying  Aq^(w)  for  the  N-element  linear  array.)  Consequently, 
using  a  set  of  taps  equally  spaced  along  a  delay  line  would  yield  a  better 
approximation  to  this  simpler  optimum  frequency  variation  (when  N^  is  con¬ 
strained  to  be  small)  than  it  would  for  the  more  complex  frequency  variation 
associated  with  the  more  complex  arrays.  Assuming  this  conjecture  to  be  valid, 
one  would  then  require  fewer  taps  per  element  for  an  N-element  linear  array 
than  for  the  more  complex  N-element  arrays  in  order  to  achieve  a  desired 
cancellation  performance.  Otherwise  said,  the  number  of  taps  required  depends 
more  on  the  spatial  complexity  of  the  array  rather  than  on  the  number  of  ele¬ 
ments  in  the  array.  The  following  results  offer  support  to  the  validity  of 
this  conjecture,  although  we  can  by  no  means  present  a  definitive  proof  of  its 
general  applicability. 

In  the  following,  we  consider  first  the  simplest  class  of  arrays  —  the  N- 
element  linear  array  —  followed  by  the  analysis  for  the  triangular  and  more 
complex  arrays. 

A.  The  N-Element  Linear  Array 

1.  Three  Elements:  The  three  element  linear  array  is  an  interesting 
extension  of  the  two-element  array  because  it  allows  for  the  presence  of  two 
interference  sources  located  arbitrarily  over  the  FOV,  yet  is  still  simple 
enough  so  that  some  analytical  results  can  be  obtained.  One  would  anticipate 
that  the  amplitude  and  phase  frequency  variation  of  the  optimum  weights  for  the 
three  element  array  becomes  more  complex  when  compared  to  that  for  the  two- 
element  array  (for  which  the  optimum  weights  have  only  phase  variation),  and 
hence  a  greater  degree  of  synthesis  capability  is  required  to  approximate  this 
variation.  Indeed,  this  turns  out  to  be  the  case,  but  we  shall  see  that  the 
cancellation  performance  of  the  array  still  depends  only  weakly  on  the  number 
of  elements. 

Consider  the  three-element,  equi-spaced  array  illustrated  in  Fig.  12,  and 
the  two-source  interference  scenario  defined  by  incidence  angles  0^  and  02* 
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k0Dxk 

In  keeping  with  notation  introduced  in  Section  III,  define  *  — 2 —  s^n 
and  2  =  ^0°^^  sin  ®2’  k=*’  or  w^ere  tl'e  ^  are  normalized  element 


If  we  choose  x^  =  -1, 


X2  =  0  and  x^  =  1  for  element 


Y2  2  reduce  to  0,  and  y1  ^  =  -Yj 


Y2,l  ~  ~y3,2- 


positions  (x^  <  1) , 
locations,  then  ^  ±  and 
This  simplifies  the  analysis  considerably.  It  is  instructive  to  determine  the 
optimum  weight  variation  defined  by  Eq.  (15)  for  this  array,  and  use  it  as  a 
basis  for  comparison  to  the  tapped  delay  line  synthesized  weight.  As  before, 
denote  these  optimum  weights  as  Aq^(w),  k  =  1,  2,  3.  Then  it  can  be  shown  that 


A0i(w) 


A03(w)  =  -e 


(42) 


a02 (w)  =  2  cos 


[ 


^2,2  ~  Y3,l^ 


"] 


Observe  that,  as  contrasted  with  the  phase  only  variations  of  the  optimum 
weighting  for  the  two-element  array,  the  A^^(w)  in  (42)  have  both  amplitude 
and  phase  dependence  on  frequency.  In  fact,  a  careful  examination  of  Eq.  (42) 
reveals  that  it  is  not  generally  possible  to  synthesize  AQk(w)  exactly 
using  a  tapped  delay  line  with  a  finite  number  of  taps  having  equal  tap-spacing. 
Consequently,  the  restriction  of  equal  tap-spacing  necessarily  yields  non¬ 
optimum  results  for  some  angles  of  arrival.  However,  due  to  its  simplicity 
and  practical  utility,  it  is  still  a  useful  constraint  to  impose  on  the  delay 
line  design.  Furthermore,  it  is  certainly  possible  that  when  the  array  perfor¬ 
mance  is  averaged  over  all  possible  sets  of  incidence  angles  allowable  over 
the  FOV,  uniform  tap-spacing  would  still  be  the  desired  choice. 

Consider  now  the  synthesis  of  Eq.  (42)  using  either  one,  two  or  three 
taps  per  element  for  the  three-element  array.  Assume  two  orthogonal  inter¬ 
ference  sources  located  over  the  FOV,  with  one  source  positioned  at  maximum 
scan  angle  0  .  A  second  orthogonal  source  will  be  orthogonal  to  the  first 
if  it  is  positioned  at  0^  satisfying  the  relationship 
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m 


(sin  0  sin  0O)  =  2ir/3  (43) 

z  m  z 

For  small  0  (as  for  a  satellite  at  geosynchronous  altitude),  we  have  0  -  Q_ 

m  m  L 

~~  2/3  X/D.  For  the  two  sources  so  positioned,  the  resultant  cancellation  as 

Ay  •  BBW  is  varied  (actually,  to  obtain  the  numerical  results,  we  fix  D/X, 
m 

0^,  0£  and  vary  FBW)  is  plotted  in  Fig.  10,  along  with  the  results  for  the  two 
element  array*.  For  lack  of  a  better  criterion,  we  have  chosen  tap-spacing  <j>g 
according  to  the  optimum  spacing  developed  for  the  two-element  array.  This 
spacing  tends  to  minimize  the  average  cancellation  when  considered  as  a  function 
of  all  angles  of  incidence  over  the  FOV.  Thus, 


=  .775  Ay  for  N  =  2 
0  m  t 

<pn  =  3875  Ay  for  N  =  3 

u  m  t 


(44) 


where  Ay  =  2ir  D/X  •  sin  0  .  We  have  also  chosen  ysWA„  >>  1  so  as  to  get 
m  m  MAX 

at  the  fundamental  nulling  limitation  of  the  array,  independent  of  the 
algorithm  implemented  in  the  processor.  Notice  the  close  similarity  between 
performance  obtained  for  the  two  and  three  element  array  configurations  as 
the  number  of  taps  is  varied.  In  Figs.  13a,  13b  and  13c,  we  examine  the 
frequency  dependence  of  the  delay  line  synthesized  weights  A^(w)  vs  the  optimum 
weights,  Aq^(w),  of  Eq.  (42)  using  Ay^  •  FBW  =  7.0.  We  utilize  the  polar  repre¬ 
sentation  of  the  weights  introduced  in  Section  II.  Since  each  weight  can  be 
multiplied  by  an  arbitrary  complex  function  of  frequency  without  changing  the 
results,  the  results  in  Fig.  13  are  normalized  to  A^(w)  for  each  respective 


*Note:  These  results,  and  all  of  the  following  simulations  assume  an  earth 
coverage  quiescent  radiation  pattern.  Any  other  steering  vector  would 
increase  or  decrease  the  cancellation  achieved  by  approximately  the  average 
quiescent  gain  to  the  interference  sources  relative  to  earth  coverage  gain. 
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In  keeping  with  notation  introduced  in  Section  III,  define  ^  <=  — 1 -  sin  0 
and  2  =  ^0Dxk^  S^n  ®2’  k=l>  or  •*»  wh®*6  th®  ^  are  normalized  clement 
positions  (x^  <  1) .  If  we  choose  =  -1,  X2  =  0  and  x^  *  1  for  element 
locations,  then  ±  and  Y£  2  reduce  to  0,  and  ^  =  “ >3  1  *  Y2  1  =  *^3  2‘ 

This  simplifies  the  analysis  considerably.  It  is  instructive  to  determine  the 
optimum  weight  variation  defined  by  Eq.  (15)  for  this  array,  and  use  it  as  a 
basis  for  comparison  to  the  tapped  delay  line  synthesized  weight.  As  before, 
denote  these  optimum  weights  as  Aq^(w),  k  =  1,  2,  3.  Then  it  can  be  shown  that 


koDxk 


AqI  (W>  A03^ 


A02(w)  =  2  cos 


("I3  2  '  ^3  1^ 

- > - 1 -  TJ 

2 


Observe  that,  as  contrasted  with  the  phase  only  variations  of  the  optimum 
weighting  for  the  two-element  array,  the  AQk(w)  in  (42)  have  both  amplitude 
and  phase  dependence  on  frequency.  In  fact,  a  careful  examination  of  Eq.  (42) 
reveals  that  it  is  not  generally  possible  to  synthesize  AQ^(w)  exactly 
using  a  tapped  delay  line  with  a  finite  number  of  taps  having  equal  tap-spacing. 
Consequently,  the  restriction  of  equal  tap-spacing  necessarily  yields  non¬ 
optimum  results  for  some  angles  of  arrival.  However,  due  to  its  simplicity 
and  practical  utility,  it  is  still  a  useful  constraint  to  impose  on  the  delay 
line  design.  Furthermore,  it  is  certainly  possible  that  when  the  array  perfor¬ 
mance  is  averaged  over  all  possible  sets  of  incidence  angles  allowable  over 
the  FOV,  uniform  tap-spacing  would  still  be  the  desired  choice. 

Consider  now  the  synthesis  of  Eq.  (42)  using  either  one,  two  or  three 
taps  per  element  for  the  three-element  array.  Assume  two  orthogonal  inter¬ 
ference  sources  located  over  the  FOV,  with  one  source  positioned  at  maximum 
scan  angle  0  .  A  second  orthogonal  source  will  be  orthogonal  to  the  first 
if  it  is  positioned  at  ©2  satisfying  the  relationship 
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Fig.  13.  Optimum  weighting  vs  delay  line  synthesized  weighting 
for  the  three-element  linear  array:  (a)  Optimum  weight  variation 
with  frequency,  (b)  Delay  line  synthesized  weight  variation  with 
frequency  using  two  taps/element. 
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Fig.  13  (Continued).  Optimum  weighting  vs  delay  line  synthesized 
weighting  for  the  three-element  linear  array:  (c)  Three  taps/ele 


\TT 


delay  line  case,  and  to  Aq^(w)  for  the  optimum  weight  variation.  This  allows 
for  comparison  on  a  one-one  basis.  Note,  as  obtained  from  Eq.  (42),  A^^(w) 
and  Aq2(w)  illustrated  in  Fig.  13a  are  constant  amplitude  weights  and  Aq^(w) 
has  a  cosine  variation  with  frequency  with  a  peak  amplitude  of  nearly  6  dB. 

The  value  of  the  Aq^(w)  at  w  =  1  is  indicated  in  the  figure.  Figs.  13b  and  13c 
illustrate  the  delay  line  synthesized  approximation  to  the  A^^(w)  using  2  and 
3  taps/element,  respectively.  At  band  center,  (w  =  1),  all  the  weights  compare 
closely.  As  w  increases  from  band  center,  the  delay  line  synthesized  weights 
begin  to  deviate  from  the  optimum  weights,  although  and  A^(w)  have 

constant  amplitude.  The  approximation  to  A^^(w)  is  worst  at  the  edge  of  band, 
and  it  is  here  that  the  improvement  realized  in  using  3  taps/element  over  2 
taps/element  can  be  seen.  This  is  particularly  true  of  A^(w)  for  N  =  2, 
where  the  approximation  to  A^^(w)  is  very  poor  at  the  edge  of  the  nulling 
band.  For  N  ■  3,  A^(w)  closely  approximates  Aq^(w).  In  order  to  illustrate 
how  these  deviations  in  A^(w)  from  aq^(w)  affect  the  cancellation  realized  as 
a  function  of  frequency,  we  illustrate  in  Figs.  14a  and  14b  the  cancellation 
vs  frequency  over  the  nulling  band  obtained  in  the  direction  of  each  of  the 
two  sources.  For  N  ■  1,  hardly  any  nulling  results,  whereas  the  improvement 
in  cancellation  realized  in  increasing  Nt  to  either  2  or  3  is  quite  significant. 
Observe  the  cancellation  is  worst  at  the  band  edges  consistent  with  the  devia¬ 
tion  of  A^(w)  from  Ap^(w)  in  these  regions.  The  average  cancellation  over  the 

nulling  band  can  be  obtained  using  Fig.  10  for  Ay^’FBW  =  7,  from  which  we 

obtain  C  - 20  dB  using  *=  2  and  C  ~  -40  dB  using  Nt  =  3.  Clearly  these 

values  are  consistent  with  the  cancellation  vs  frequency  results  illustrated 
in  Fig.  14. 

2.  The  Six-Element  Linear  Array:  In  order  to  obtain  a  feel  for  the 
dependence  of  the  number  of  taps  on  the  number  of  elements,  consider  the 
six-element  array  illustrated  in  Fig.  11a.  For  this  array,  five  orthogonally 
located  interference  sources  can  be  positioned  nearly  uniformly  over  the  FOV, 
similar  to  the  positioning  illustrated  shown  in  the  insert  of  Fig.  1.  As 
before,  we  consider  the  cases  N  ■  1,  2  or  3  using  *  .775  A  ^  for  Nfc  =  2 

and  4> _  =  0.3875  Ay  for  =  3.  The  resultant  cancellation  performance  when 

u  m  t 
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plotted  vs  Af^'FBW  is  illustrated  in  Fig.  10,  for  comparisor  with  the  two  and 

three  element  array  results.  The  general  similarity  of  the  results  for  all 

these  arrays  when  -  2  and  »  3  indicates  that  the  performance  achieved 

is  only  weakly  dependent  on  the  number  of  elements.  Regions  where  one,  two 

or  three  taps  should  be  used  to  obtain  a  desired  cancellation  performance  are 

readily  identifiable.  Note  also  the  general  applicability  of  the  cancellation 

2 

estimates  obtained  from  the  two-element  array  which  predict  C  «  (Ay *FBW)  for 

6  10  ™ 

N  ■  1,  C  «  (Ay^'FBW)  for  Nfc  =  2  and  C  «  (Ay^'FBW)  for  N  *  3.  Hence,  we 

conclude  that  for  the  linear,  equi-spaced  array  the  key  parameters  governing 

the  array  performance  are  D/A ,  sin  6  and  the  fractional  bandwidth.  To  first 

zn 

order  these  parameters  can  be  used  to  estimate  the  number  of  taps  required  in 
order  to  achieve  a  desired  wideband  cancellation  performance. 

B.  Triangular  Arrays 

The  triangular  arrays  illustrated  in  Fig.  lib  represent  the  simplest 
point  of  departure  in  considering  the  transition  from  linear  arrays  to  planar 
arrays.  In  treating  even  simple  planar  arrays,  the  analysis  becomes  complicated 
due  to  the  necessary  inclusion  of  a  two-dimensional  FOV.  However,  the  general 
concept  of  orthogonally  positioned  interference  sources  and  their  complete  use 
of  an  array  degree  of  freedom  remains  valid.  Consequently,  the  general  phil¬ 
osophy  of  positioning  N-l  orthogonal  (or  nearly  so)  sources  and  examining  the 
cancellation  performance  of  the  array  relative  to  these  sources  as  the  nulling 
bandwidth  increases  can  still  be  used  to  evaluate  the  array  performance.  Since 
aperture  frequency  dispersion  is  worse  for  sources  located  toward  the  edge  of 

the  FOV  than  for  those  near  the  center  (see  Fig.  1,  cancellation  without 

2 

frequency  compensation  degrades  as  (A^FBW)  ),  we  consistently  position  N-l 
sources  at  these  locations  whenever  possible.  It  should  be  emphasized,  however, 
that  when  0/A  is  large  enough  so  that  grating  lobes  occur  over  the  FOV,  the  set 
of  N-l  orthogonal  source  locations  are  not  unique  and  cancellation  performance 
is  dependent  on  the  set  of  source  locations  chosen.  Consequently,  locations 
chosen  in  the  following  analysis  yield  results  which  are  only  representative 
of  the  bandwidth  characteristics  of  the  array,  and  are  not  necessarily  worst 
case. 
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With  these  restrictions  in  mind,  consider  the  performance  of  the  4-element 
triangular  array  Illustrated  in  Fig.  lib,  subjected  to  three  interference 
sources  orthogonally  located  over  the  FOV.  The  cancellation  vs  Ay^'FBW  is 
illustrated  in  Fig.  15  for  two  different  source  location  scenarios,  illustrated 
in  the  insert  in  the  figure,  using  values  of  N  ■  1,  2  and  3.  The  performance 
of  the  array  vs  the  two  scenarios  is  somewhat  different  for  each  case.  However, 
even  the  poorest  performance  yields  significant  cancellation  over  the  nulling 
band  as  the  number  of  taps  increases.  Perhaps  most  significant  is  that  using 
two  taps  per  element  yields  only  about  20-30  dB  improvement  over  frequency 
independent  weighting,  whereas  using  3  taps  per  element  results  in  better  than 
50  dB  (theoretical)  improvement  over  that  obtained  using  frequency  independent 
weighting.  Furthermore,  the  large  cancellation  levels  realized  tend  to  justify 
the  choice  of  4> Q  used,  which  was  extrapolated  without  justification  from  the 
two-element  array  analysis. 

In  order  to  investigate  the  dependence  of  cancellation  on  the  number  of 
elements,  we  fix  the  array  symmetry  (i.e.,  3  radial  arms  specifying  the  element 
locations)  and  increase  the  number  of  elements  from  four  to  seven.  The  cancel¬ 
lation  vs  Ay^-FBW  for  this  array  ia  illustrated  in  Fig,  16  for  the  6-orthogonal 
source  scenario  Indicated  in  the  figure.  Clearly,  the  cancellation  realized 
as  a  function  of  A^'FBW  is  not  markedly  different  from  that  for  the  four- 
element  triangular  array,  again  indicating  only  a  weak  dependence  of  the  delay 
line  synthesis  of  A^^(w)  on  the  number  of  elements.  In  fact,  recognizing  that 
the  triangular  array  is  comprised  of  three  linear  arrays  rotated  120°,  it 
should  not  be  surprising  that  the  results  are  not  markedly  different  from  the 
linear  array.  This  supports  the  conjecture  that  the  number  of  taps  required 
in  order  to  achieve  very  high  cancellation  levels  is  most  strongly  dependent 
on  the  angular  dispersiveness  of  the  array,  i.e.,  required  is  proportional 
to  the  number  of  radial  arms  characterizing  the  array. 

C.  More  Complex  Arrays 

Consider  now  the  class  of  arrays  illustrated  in  Fig.  11c,  i.e.,  those 
requiring  more  than  three  angularly  placed  radial  arms  to  specify  their  element 
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CANCELLATION  (dB) 


Fig.  16.  Cancellation  (dB)  vs  Ay  •  FBW  for  the  seven-element 
simple-triangular  array  using  a  s¥x-source  orthogonal  scenario  as 
a  function  of  the  number  of  taps. 


locations.  In  order  to  illustrate  the  effect  of  radial  element  positioning  on 
the  cancellation  performance,  consider  a  simple  modification  to  the  7-element 
triangular  array  treated  above.  If  the  inner  triad  of  array  elements  of  the 
7-element  triangular  array  of  Fig.  lib  is  rotated  60°,  the  "rotated-double 
triangle"  array  illustrated  in  Fig.  11c  results.  This  simple  rotation  of 
elements  results  in  an  array  configuration  characterized  by  6  equi-spaced 
radial  arms.  The  resultant  cancellation  vs  A-^*FBW  performance  relative  to 
nulling  a  six- interference  source  scenario  is  illustrated  in  Fig.  17  using  a 
one,  two  or  three  tapped  delay  line.  The  six-source  scenario  used  to  evaluate 
the  array  performance  is  a  nearly-mutually  orthogonal  scenario  (8  dB  spread 
between  the  first  six  eigenvalues),  chosen  because  a  set  of  six  mutually 
orthogonal  sources  for  this  array  configuration  does  not  appear  to  exist. 
Comparing  Fig.  17  to  Fig.  16,  we  conclude  that  the  simple  notational  change 
from  three  to  six  radial  elements  drastically  alters  the  array  nulling  band¬ 
width  performance.  This  is  because  the  more  spatially  dispersive  the  array, 
the  more  difficult  it  is  for  a  fixed  number  of  taps  on  a  delay  line  to  synthe¬ 
size  the  optimum  weight  variation  with  frequency.  To  see  this,  we  illustrate 
in  Figs.  18  and  19  the  optimum  synthesized  tapped  delay  line  weights  variation 
with  frequency  for  the  scenario  considered  using  the  simple  double  triangle 
configuration  of  Fig.  lib  and  the  rotated  double  triangle  configuration  of  Fig. 

11c,  respectively,  and  a  value  of  Ay  *FBW  =  3.43.  For  each  figure,  the  polar 

m 

plot  labeled  "a"  illustrates  the  optimum  variation  with  frequency  and  the  plots 
labeled  "b"  or  "c"  illustrate  the  tapped  delay  line  synthesized  weights  using 
two  and  three  taps/element,  respectively.  Note  the  rather  simple  amplitude  and 
phase  variation  of  the  optimum  weights  with  frequency  for  the  simple  double 
triangle  array,  as  compared  to  the  rather  complex  variation  obtained  for  the 
rotated  double  triangle.  For  the  simple  double  triangle  array,  the  resultant 
optimum  and  synthesized  weights  lie  on  the  real  axis  for  the  scenario  chosen. 
Because  of  this,  it  is  difficult  to  assess  from  the  polar  plot  quantitative 
differences  between  the  optimum  weights  and  the  results  for  Nt  =  2  and  N  ■*  3, 
as  the  data  is  compressed  onto  a  line.  However,  a  more  detailed  examination  of 
this  weight  frequency  variation  shows  a  nearly  monotonic  variation  with 
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Fig.  18.  Optimum  weighting  va  delay  line  synthesized  weighting 
for  the  seven-element  simple-double  triangular  array  using  Ay 
FBW  =  3.43:  (a)  Optimum  weight  variation  with  frequency,  (b)  mDelay 

line  synthesized  weight  variation  with  frequency  for  *  2  and  3. 
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Fig.  19.  Optimum  weighting  vs  delay  line  synthesized  weighting 
for  the  seven-element  rotated  triangular  array  for  Ay  •  FBW  *  3.43 
(a)  Optimum  weight  variation  with  frequency,  (b)  Delay  line  synthe 
sized  weight  variation  with  frequency  using  N  =  2. 


frequency  results  for  the  simple  double  triangle  array,  as  opposed  to  that 
corresponding  to  the  rotated  double  triangle,  for  which  some  rather  abrupt 
Inflection  points  occur  for  the  phase  variation  with  frequency.  As  can  be  seen 
in  Figs.  19b  and  c,  using  only  two  or  three  taps  is  not  sufficient  to  precisely 
synthesize  the  rapid  phase  variation  of  the  optimum  weights  of  Fig.  19a, 
although  the  approximation  to  the  optimum  weights  clearly  improves  as  the 
number  of  taps  increases.  One  might  also  compare  the  optimum  weight  variation 
with  frequency  for  these  two  triangular  array  configurations  to  that  of  Fig.  13 
for  the  three-element  linear  array,  where  the  smooth  behavior  of  the  phase 
variation  is  again  evident. 

In  order  to  pursue  this  point  further,  consider  another  7-element  array 
configuration  characterized  by  six  radial  arms  —  the  hexagonal  array  illus¬ 
trated  in  Fig.  11c.  Six  orthogonal  source  locations  for  this  array  can  readily 
be  determined  and  are  illustrated  in  the  insert  of  Fig.  20.  The  cancellation 
performance  vs  the  parameter  Ay^'FBW  is  also  plotted  in  this  figure  for  a  delay 
line  employing  either  one,  two  or  three  taps/element.  The  tap  spacing  is,  as 

before,  chosen  to  be  4>^,  *»  .775  Ay  for  =  2  and  .3875  Ay  for  N  =  3.  We  note 

u  m  t  m 

the  relatively  narrowband  nulling  performance  of  this  array  is  slightly  worse 
than  for  the  7-element  rotated  double  triangle  array.  As  one  might  anticipate 

from  the  preceding  results,  the  optimum  weight  variation  with  freqijeney  also 

/ 

becomes  rapidly  varying  as  Ay  •  FBW  increases.  This  is  illustrated  in  Fig.  21, 

m 

where  we  plot  this  weight  variation  vs  frequency  for  A-^FBW  =  2.19.  The 
optimum  weight  variation  with  frequency  is  plotted  in  Fig.  21a,  along  with  the 
synthesized  weight  variation  with  frequency  in  Figs.  21b  and  c,  respectively. 
Clearly,  the  same  general  comments  pertaining  to  the  rotated  double  triangle 
relative  to  Fig.  19  can  be  made  for  the  hexagonal  array.  The  inflection 
points  present  in  the  optimum  weights  phase  variation  occurs  both  near  band- 
center  and  at  the  band  edges. 

For  completeness,  we  briefly  consider  the  two  remaining  array  configura¬ 
tions  of  Fig.  11c:  The  square  array  and  the  pentagonal  array,  characterized  by 
four  and  five  radial  arms,  respectively.  The  cancellation  performance  for 
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Fig.  20.  Cancellation  (dB)  vs  Ay  •  FBW  for  the  seven-element 
hexagonal  array  using  a  six-source  orthogonal  scenario  as  a  function 
of  the  number  of  taps. 


♦ 
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Fig.  21.  Optimum  weighting  vs  delay  line  synthesized  weighting  for 
the  seven-element  hexagonal  array  for  Ay  *  FBW  »  2.19:  (a)  Optimum 
weight  variation  with  frequency,  (b)  Delay  line  synthesized  weight 
variation  with  frequency  using  Nt  ■  2  . 
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these  two  arrays  vs  A^'FBW  is  illustrated  in  Figs.  22  and  23,  respectively. 

The  nearly  orthogonal  source  locations  used  in  the  simulation  are  illustrated 
in  the  insert  of  each  figure.  For  comparison,  we  have  presented  results  for 
two  different  orthogonal  source  scenarios  for  the  5-element  square  array,  and 
included  results  using  4  taps/element  for  both  arrays.  When  *  4,  we  choose 
4>q  *  A  “tjjj/4.  The  results  indicate  performance  similar  to  the  other  array  types 
of  Fig.  11c.  Observe  that  notable  improvement  is  obtained  when  4  taps/element 
are  used,  and  the  results  become  comparable  to  the  more  simple  linear  and  tri¬ 
angular  arrays  using  two  delay  lines/element. 

Finally,  in  order  to  bring  together  the  results  for  all  of  the  array  types 
in  Fig.  11,  we  examine  the  composite  performance  of  the  arrays  obtained  using 
a  fixed  number  of  taps /element.  To  accomplish  this,  we  Illustrate  in  Fig.  24a 
the  collective  performance  of  each  array  type  vs  Ay^FBW  using  frequency 
independent  weights.  Also  sketched  in  the  figure  is  the  simple  "rule  of  thumb" 
discussed  in  Section  1  (Fig.  1)  used  to  assess  the  need  for  employing  tapped 
delay  lines.  Note  this  rule  of  thumb  curve  serves  as  a  useful  upper  bound  on 
the  composite  set  of  cancellation  results  for  the  general  class  of  arrays 

2 

considered.  Observe  that  for  all  of  the  array  configurations,  C  <*  (Ay^FBW) 

--  i.e.,  for  frequency  independent  weighting,  array  bandwidth  performance 

degrades  as  the  square  of  the  parameter  2it*^  sin  0  *FBW.  This  is  a  useful 

a  in 

relationship  to  keep  in  mind  in  assessing  array  performance. 

Fig.  24b  illustrates  the  collective  improvement  of  the  9  array  types 
considered  using  >  2.  Observe  a  much  larger  spread  in  cancellation  perfor¬ 
mance  is  obtained,  indicating  the  Improvement  realized  is  somewhat  sensitive  to 
array  configuration.  However,  the  spread  is  narrow  enough  to  enable  one  to 
estimate  useful  bounds  on  the  anticipated  improvement  realized  using  *  2  as 

a  function  of  A^'FBW.  These  bounds  can  be  tabulated  as  a  function  of  the 
desired  cancellation  level  for  a  given  number  of  taps.  (Recall,  all  the  above 
results  were  obtained  assuming  an  earth  coverage  quiescent  steering  vector. 

Any  other  steering  vector  would  increase  or  decrease  the  cancellation  achieved 
by  approximately  the  average  quiescent  gain  to  the  interference  sources  rela¬ 
tive  to  earth  coverage.)  The  result  is  illustrated  in  Table  I.  Fig.  24c 


Fig.  22.  Cancellation  (dB)  vs  Ay  •  FBW  for  the  five— element  square 
array  using  two-different  four-source,  nearly  orthogonal  scenarios 
as  a  function  of  the  number  of  taps. 
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Fig.  2A.  Composite  performance  of  the  nine  array  configuration 
illustrated  in  Fig.  11:  (a)  Using  frequency  independent  weightl 
(b)  Using  delay  line  weighting  with  two  taps/element. 
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Fig.  24  (Continued) .  Composite  performance  of  the  nine  array  con¬ 
figurations  illustrated  in  Fig.  11:  (c)  Using  delay  line  weighting 

with  three  taps/element. 


TABLE  1 


COMPOSITE  CANCELLATION  BOUNDS  vs  N  AND  X  -  Ay  -FBW  FOR  THE  NINE  ARRAY 

t  m 

CONFIGURATIONS  OF  FIG.  11  ASSUMING  AN  EARTH  COVERAGE  QUIESCENT 

RADIATION  PATTERN 
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illustrates  the  composite  cancellation  performance  using  N  ■  3.  Observe  that 
those  classes  of  arrays  categorized  in  Fig.  11  yield  much  different  performance 
characteristics,  offering  support  to  the  conjecture  relating  cancellation 
dependence  on  the  number  of  taps  to  the  spatial  distributions  of  elements  over 
the  aperture.  Results  for  *  3  tabulated  in  Table  1  use  the  simulation  data 
for  the  more  complex  arrays,  i.e.,  the  upper  bound  of  Fig.  24c,  which  serves  as 
a  worst  case  estimate.  Although  the  bounds  are  somewhat  loose,  they  should 
prove  useful  in  estimating  a  given  array  performance  based  on  the  given  D/A, 
field  of  view,  and  operating  bandwidth. 
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V.  DISCUSSIONS  AND  CONCLUSIONS 

The  results  of  Section  IV  indicate  the  cancellation  performance  of  various 
types  of  planar  arrays.  It  was  determined  that  the  cancellation  potential  of 
a  particular  array  could  be  considered  a  function  of  the  parameter  Aym*FBW. 

Some  physical  insight  into  the  interpretation  of  this  parameter  can  be  obtained 
in  a  qualitative  manner  by  consideration  of  the  time  delays  associated  with  the 
delay  lines  and  the  array.  The  maximum  delay  between  the  times  of  arrival  of 
a  plane  wave  at  any  pair  of  elements  in  an  array  of  diameter  D  is  given  by 

D  sinG 


where  c  is  the  velocity  of  light.  A  perfect  broadband  weight  would  be  capable 

of  compensating  for  any  time-delay  from  0  to  (see  Eq.  (18)).  The  choice 

of  a  total  delay  line  length  which  gives  a  time  delay  of  775  x.,.„  permits 

‘  MAX 

one,  by  appropriate  weighting  of  the  taps  along  the  line,  to  approximate 
optimum  weighting  over  this  range  relative  to  an  arbitrary  delay  over  the 


range  0  <  t  <  t 


Simulations  have  shown  the  nulling  performance  to  be 


not  too  sensitive  to  the  value  of  total  delay  line  length.  For  the  case  of 

*  ^  ^UAY 

3  taps,  we  have  somewhat  arbitrarily  chosen  the  tap  spacing  to  be  _ _ ir** 


=  .3875  For  a  larger  number  of  taps,  however,  the  appropriate  spacing 

would  appear  to  be  TMAX/Nt»  Nt  being  the  number  of  taps.  The  choice  of 
uniform  tap  spacing,  while  somewhat  arbitrary,  appears  reasonable,  given  a 
fixed  maximum  length  for  the  delay  line.  Such  a  uniform  distribution  of  taps 
should  permit  most  close  approximation  to  any  desired  delay  between  0  and 

TMAX' 

The  number  of  taps  required  to  achieve  a  certain  cancellation  ratio 
depends  on  and  on  the  fractional  bandwidth.  This  dependence  can  be 

clarified  by  introducing  the  concept  of  a  characteristic  time,  Tq,  defined 
as  the  inverse  of  the  actual  bandwidth  in  hertz;  i.e.. 


/**•«  '* 
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1 

T  *  - 

0  BW 


(46) 


or,  in  terms  of  the  fractional  bandwidth  FBW, 


T 


0 


1 

c  •  FBW 


(47) 


The  parameter  iQ  may  be  thought  of,  for  example,  as  the  width  of  the  envelope 
of  an  RF  pulse,  the  major  frequency  components  of  which  lie  within  the  spectrum 
BW.  It  is  clear  that  in  order  for  cancellation  of  such  a  pulsed  waveform 
to  occur,  the  pulse  must  appear  simultaneously  at  the  output  of  taps  on  more 
than  one  array  element.  For  a  two-element  array,  for  example,  the  number 
of  taps  along  a  delay  line  of  total  length  t  must  be  much  greater  than 

tmax 

”‘T  '  in  order  to  be  able  to  find  such  a  coincidence  of  pulses  at  the  array 

output  for  any  angle  of  arrival  up  to  0^. 

(Using  (45)  and  (47)),  the  parameter  TMA^/T g  can  expressed  as 


2ir 


MAX 


=A  y 


m 


FBW 


(48) 


Eq.  (48),  when  interpreted  relative  to  tj^ax/to»  °fferfl  insight  into  the 

physical  significance  of  the  parameter  Ay^  •  FBW  used  in  plotting  many  of  the 

curves  in  the  report.  For  example,  if  only  a  single  tap  is  used,  it  is 

reasonable  to  require  <  •!»  or  A •  FBW  <  .63,  and  the  curves  in 

the  report  show  poor  cancellation  for  Ay  •  FBW  >.63. 

tn 

Other  points  that  the  authors  wish  to  reiterate  are  that  our  analysis 
always  assumes  a  maximally  stressed  array,  i.e.,  an  N-element  array  is  always 
evaluated  relative  to  N-l  orthogonal  sources  distributed  over  the  FOV.  This 
assumes  that  the  array  dedicates  only  a  single  spatial  degree  of  freedom  per 
source  to  be  nulled. 


A  different  approach  to  obtaining  the  capability  of  nulling  N-l  sources 
over  a  wide  band  might  be  to  use  M  >>  N  elements  distributed  over  the  same 
diameter  D,  but  no  tapped  delay  lines,  so  that  many  degrees  of  freedom  could 
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be  devoted  to  wideband  nulling  of  each  source.  The  disadvantage  of  this 

approach  is  that  the  nulls  obtained  are  broader  in  angle,  so  that  the  nulling 

/ 

resolution  is  degraded,  an  effect  that  does  not  occur  with  the  tapped  delay¬ 
line  approach.  It  also  requires  the  physical  deployment  of  many  more  elements 
than  would  be  necessary  with  tapped  delay  lines.  The  tapped-delay  line  permits 
one  to  approximate  the  ideal  narrowband  weight  at  each  frequency,  and  thus  to 
maintain  a  high-resolution  null  (consistent  with  the  resolving  power  of  the 
antenna  aperture)  in  a  desired  direction. 

We  have  also  shown  that  the  number  of  taps  required  per  array  element 
depends  not  only  on  the  parameter  A  •  FBW,  but  also  on  the  geometry  or  shape 
of  the  array,  more  taps  being  required  for  more  complex  shapes.  The  number 
of  taps  is  only  weakly  dependent  on  the  number  of  elements  that  are  placed 
inside  the  diameter  D. 

The  analysis  has  been  carried  out  using  the  Applebaum-Howells  algorithm 
to  obtain  optimum  tap  weighting,  but  care  has  been  taken  to  make  the  results 
of  Section  IV  independent  of  the  particular  algorithm  used.  This  has  been 
done  by  choosing  a  large  enough  value  of  y,  the  adaptive  processor  loop  gain, 
and  by  choosing  an  appropriate  quiescent  weight  vector,  one  which  is  guaranteed 
not  to  be  orthogonal  to  the  desired  nulling  weight  vector. 

In  Table  I  in  the  last  section  of  the  report,  we  attempted  to  summarize 
the  bandwidths  achievable  for  1,  2  and  3  taps  at  different  cancellation  levels. 
The  table  lists  a  range  of  values  because  of  the  above  mentioned  effect  of 
dependence  of  the  results  on  the  shape  of  the  array.  This  table  can  be  used 
as  a  preliminary  design  guide,  before  more  detailed  analysis  is  carried  out. 


APPENDIX  A 


AN  UPPER  BOUND  ON  THE  CANCELLATION  ACHIEVABLE 
USING  FREQUENCY  INDEPENDENT  WEIGHTING 

In  this  appendix,  we  develop  the  results  leading  up  to  Fig.  1  of  the  text, 
which  serves  as  a  qualitative  estimate  of  the  bandwidth  potential  of  the  N- 
element  array  antenna  using  frequency  independent  weights.  Assume  the  presence 
of  N-l  orthogonally  located  sources.  Then  Eq.  (4)  characterizing  R  can  be 
expanded  in  the  form 


S  -  So  +  AR  =  Si.o  +  •••  +  Bn-1,0  +  A£i  +  —  +  aSn-i 


CAl) 


where  R.  _  denotes  the  zero-bandwidth  correlation  matrix  of  the  jth  source  and 

AR.  the  perturbation  to  R  _  due  to  non-zero  bandwidth.  The  Nth  eigenvalue 
— 3  j  » u 

sN  is  given  by 


SN  =  MaX  - 1 - 

e  e  •  e 


R*£  £  •  AR  •  e 

Z -  =  Max  - ^ - 

e  e  •  e 


(A2) 


where  the  maximization  in  (A2)  is  performed  over  all  e  orthogonal  to  the  N-l 
eigenvectors  of  Rq.  Using  a  simple  inequality,  it  follows  that  sN  can  be 
bounded  by 


N-l 


SN<  2 

J-l 


(A3) 


where  no  restriction  is  imposed  on  e.  Hence  the  maximum  eigenvalue  of  each 

AR.  can.be  estimated,  and  s„  bounded  by  their  sum.  Consider  AR,  and  define 
=j  kQD  N  =1 

^  sinO^  (x^cos4>^  +  y^sin^).  Then  using  Eq.  (19)  of  the  text,  it 

follows  that 


^0  2  2  ^  \  ' 
14  FB"  <\,1  -  >,.!>  e 


Vi' 


(A4) 
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As^,  corresponding  to  the  maximum  eigenvalue  of  AR^,  can  be  estimated  by  a 
simple  bound  on  the  elements  of  AR^.  From  Reference  4,  it  can  be  shown  that 

E  2  N  N 

AS1<  “§4  FBw2  (T  sinQi)2  ^  Z  2  "  *q)2  ekeq* 

k=l  q=l 


(A5) 


<  E. 


N  N 


~  •  FBW2  (—  sin01)2MAX  2  2  (\  -  xq>' 

k=l  q=l 


If  we  assume  an  N-element,  equi-spaced  linear  array  geometry,  then  the  maximum 
sum  in  (A5)  can  be  determined.  We  obtain 


MAV  V  V  /  ^2  _  2  N(2N-1) 

MAX  ^  2j<-\  ~  xq^  3  N-l 

k  q 


(A6) 


Hence 


0  2 
Asi<  FBW 


,ttA  .  „  ,2  2  N(2N-1) 
C-  sinep  3  -^33— 


(A7) 


Then  it  follows  from  (A3)  that  s^  <  As^  +  ...  +  AsN_^,  so  that 
E  2 

.  0  1  N(2N-1) 


N  24  6 


(2tt  •  ~  •  FBW)  •  [sin  ©3  +  ...  +  sin  0^] 


(A8) 


Assume  an  earth  coverage  quiescent  radiation  pattern.  The  the  quiescent 

interference-thermal  noise  ratio  is  simple  (N-l)G  R.  if  all  N-l  interference 

e  j 

sources  are  of  equal  power  level.  The  adapted  interference-thermal  noise 
ratio  is  simply  s^.  Hence,  the  cancellation  C  is  given  by 

C  s.,/ (N-l)G  R.  (A9) 

N  e  j 
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Using  Gq.  (A8)  in  (A9)  results  in  an  upper  bound  on  the  achievable  cancellation 


C  < 


_1 _ 

16ir2 


N  -v  ■<2N:1)  (2tt 

(N-l) 2 


[sin20^  +  ...  +  sin2©N 

(A10) 


where  the  ©^ ,  j  =  1,2,  N-l,  denote  the  positions  of  N-l  orthogonally 

positioned  sources.  Although  the  bound  in  Eq.  (A10)  is  not  generally  tight, 
it  serves  as  a  worst  case  estimate  of  the  bandwidth  potential  of  the  array. 

If  the  array  diameter  D/X  is  large,  so  that  grating  lobes  appear  over  the 
antenna  FOV,  then  there  are  many  angles  0.  which  satisfy  the  orthogonality 
condition  defined  by  Eq.  (2).  Since  aperture  frequency  dispersion  is  worse 
for  sources  positioned  toward  the  edge  of  the  FOV,  we  can  obtain  a  good 
estimate  of  the  worst  case  cancellation  by  choosing  the  O.'s  uniformly  distri¬ 
buted  over  the  FOV,  with  one  at  the  edge  and  the  others  at  orthogonally 
located  positions.  In  this  case,  Eq.  (A10)  can  be  expressed  in  the  form 


^  1 _  N(2N-1) 

^  o  o 

16-rr  (N-l; 


<2tt  • 


r-  sin© 

X  m 


-  sin  0 

FBW;  [1  +  - ^  +  ••• 

sin  0 

m 


sin2© 


N-2 


sin2© 


] 


(All) 


If  no  grating  lobes  occur  over  the  FOV,  then  the  orthogonal  positions  0^,  ..., 
©N_^  are  uniquely  specified;  in  this  case  Eq.  (All)  will  be  weakly  dependent 
on  N.  Fig.  1  illustrates  the  bound  of  Eq.  (All)  for  various  values  of  N  when 
the  N-l  orthogonal  source  locations  are  assumed  uniformly  distributed  over  the 
FOV.  The  tightness  of  the  bound  in  (All)  can  be  seen  by  comparison  with  simu¬ 
lations  of  the  nine  array  configurations  summarized  in  Fig.  24a. 


APPENDIX  B 


EIGENVALUE  AND  EIGENVECTOR  ANALYSIS  FOR  THE  TWO- ELEMENT  ARRAY 
HAVING  TWO  TAPS/ELEMENT 

In  this  Appendix,  we  outline  the  perturbation  analysis  leading  up  to 
Eqs.  (23)  and  (24)  characterizing  the  two-element  array  having  two  taps / 
element.  The  si  and  e^,  i  =  1,  ...»  4  defined  by  Eq.  (22)  of  the  text  can 
be  estimated  from  s^  q  and  e^  q,  which  are  the  eigenvalues  and  eigenvectors 
of  Rq,  respectively.  However,  determining  q»  o  ant*  ^4  0  comP^^cate<^ 
by  the  fact  that  the  eigenvalue  of  Rq  are  degenerate;  i.e.,  q  “  2Eq  and 
S2  0  *  s3  0  °  fl4  0  =  O'  Consequently,  the  corresponding  eigenvectors  of  these 
degenerate  eigenvalues  are  not  unique.  To  circumvent  this  problem,  we  elim¬ 
inate  the  ambiguity  in  the  following  manner.  Define  u^  q,  u^  q,  u^  q  and  H*,o 
to  be  a  given  orthonormal  set  of  eigenvectors  of  Rq.  The  desired  set  of  eigen¬ 
vectors  (e^  q)  can  be  expressed  in  terms  of  the  u^  q  with  the  relating 
constants  evaluated  using  Eq.  (22a)  of  the  text.  Clearly,  we  must  have  u^  q  = 

— 1  0*  —2  0  c^osen  so  t^at  s2  =  0  *  AR  *  — 2  0^*  where  —2  0  is 
normalized  to  unity  and  e2  o  *  — 1  0  =  Given  £2  O’  — 3  0  and  ^4  0  Can  be 
determined  using  a  similar  procedure.  Once  e^  Q  and  e^  Q  are  determined,  then 
the  e^  can  be  estimated  from  Eq.  (22b).  To  this  end,  the  eigenvector  e 2  q  is 
expanded  in  the  form 


“2,0 


=  u 


2,0 


+k3^3,0+k4^4,0 


(Bl) 


and  the  complex  constants  k^  and  k^  are  chosen  so  as  to  maximize  the  expression 
et  n  •  AR  •  e7  n 

MAX  ^  (B2) 

k3*k4  —2,0  '  -2,0 


£2  q  determined  in  this  manner  then  yields  directly  the  estimate  s^  for  the 

second  eigenvalue  of  R  given  by  Eq.  (27a).  Once  q  is  determined,  e^  q  and 

e.  .  can  be  determined  from  the  linear  combination 
-4,0 
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I 


—3,0  “  d2  —2,0  +  d3  -3,0  +  d4  ^4,0 

(B3) 

^4,0  *  c2  —2,0  +  c3  —3,0  +  c4  %0 

where  c^,  c2»  and  d^,  d2»  are  the  complex  constants,  c^,  C2»  c^,  d^, 
d2  and  d^  are  not  Independent  but  must  satisfy  the  orthonormal  conditions 


-2,0  *  -3,0  *  -3,0  '  ^4,0  =  ^4,0  *  £2,0  =  ° 


(B4) 


t  t  + 

P  •  p  ■  p  •  p  z  p  •  p  x  1 

-2,0  -2,0  -3,0  -3,0  ^4,0  ^4,0 

Eq.  (B4)  leads  to  the  set  of  linear  constraints 

c2  +  k*  c3  +  k*  c4  =  0 

d2  +  k*  d3  +  k*  d4  -  0  (B5) 

•ff  *Hf 

c2  d2  +  c3  d3  +  c4  d4  "  0 

This  set  of  equations  can  be  simplified  by  arbitrarily  choosing  d2  =  0,  so 

that  e3  q  is  dependent  only  on  u3  q  and  q  (that  this  is  possible  was 

determined  in  retrospect  by  considering  d2  ^  0,  and  is  due  to  our  choice  of 

the  u^  q) .  After  some  manipulation,  the  expressions  (B3)  for  e3  ^  and  e^  ^ 

can  be  written  in  terms  of  u,  u,  .  and  u,  in  the  form 

— 2 , 0  — 3 , 0  —4,0 


"3,0 


n  r  jy2  jYi  l 

T  [e  ^4,0  “3,0  J 

T  jY2  jYl  *o  1 

^»°  y/?  &  [S  ^.0  +  e  -3,0  A"y  “2 , 0 J 


VT 

"J  Yo 


(B6) 


(B7) 


rlw  I 
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/  2  2 

where  A  =V4>g  My  .  Eqs.  (B6)  and  (B7)  now  form  an  orthonormal  set  of  eigen¬ 
vectors  from  which  s^,  s^  can  be  determined.  It  remains  to  choose  the  basis 

set  {u.  „} .  For  our  choice  of  basis  functions,  we  choose 

—1,0 

t  i  -jYj_  -JY2  j^g  -jYx  J*0  "jY2 

-1,0  =2^e  »  e  ,e  ®  *  e  6  ^ 

t  e'jTl  „  +»•<  +«0 

u,  „  =  — 5 —  [1,  -e  ,  e  ,  -e  e  ] 


t  1  Jy0 

u I  n  =  -  [1,  0,  -e  U,  0] 

3’°  \/T 

+  1  rn  i  n  ~J*0. 

u,  n  s  -  [0,  1,  0,  -e  ] 


Then  Eq.  (B6)  for  e^  ^  used  In  (27a)  yields  and  Eq.  (B7)  for  e^  ^  gives  the 
estimate  of  s^.  Although  this  is  not  immediately  clear,  it  can  be  shown  by 
direct  maximization  of  AR  of  a  linear  combination  e  =  e_  .  =  cie,  where  a 
is  the  complex  maximization  parameter.  Having  determined  %  q*  —2  0’  —3  0 
and  e^  ^  and  s^,  S2,  s^  and  s^,  Eq.  (27b)  can  be  used  to  determine  e^,  £g 
_e^  and  e^.  However,  due  to  the  choice  of  q,  eg  q»  *.3  q  and  q»  t'ie  sum 
in  (22b)  is  zero  since  terms  i“j  are  omitted.  That  is,  e.  +  -  •  AR  •  e.  n=  0 

J  9  ^  £  *■»  V 

when  j  ^  i,  which  is  valid  for  all  terms  of  AR  up  to  order  FBW  .  Consequently, 
we  have 


-1*  -1,0’  -2  *  —2,0’  -3  =»  —3,0’  ^4  *^4,0 


where  the  e^  q  are  determined  as  above  by  substituting  (B8)  into  (B7).  The 
results  lead  to  Eq.  (24)  of  the  text. 
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